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Abstract. We combine Bethc Ansatz and field theory methods to study the 
longitudinal dynamical structure factor S zz (q,ui) for the anisotropic spin-1/2 chain 
in the gaplcss regime. Using bosonization, we derive a low energy effective model, 
including the leading irrelevant operators (band curvature terms) which account for 
boson decay processes. The coupling constants of the effective model for finite 
anisotropy and finite magnetic field are determined exactly by comparison with 
corrections to thermodynamic quantities calculated by Bethe Ansatz. We show that a 
good approximation for the shape of the on-shcll peak of S zz (q,to) in the interacting 
case is obtained by rescaling the result for free fermions by certain coefficients extracted 
from the effective Hamiltonian. In particular, the width of the on-shcll peak is argued 
to scale like Sui q ~ q 2 and this prediction is shown to agree with the width of the 
two-particle continuum at finite fields calculated from the Bethe Ansatz equations. 
An exception to the q 2 scaling is found at finite field and large anisotropy parameter 
(near the isotropic point). We also present the calculation of the high-frequency tail 
of S zz (<?, u>) in the region Suj q <C u> — vq <C J using finite-order perturbation theory in 
the band curvature terms. Both the width of the on-shcll peak and the high-frequency 
tail are compared with S zz (q, u>) calculated by Bethe Ansatz for finite chains using 
determinant expressions for the form factors and excellent agreement is obtained. 
Finally, the accuracy of the form factors is checked against the exact first moment 
sum rule and the static structure factor calculated by Density Matrix Renormalization 
Group (DMRG). 
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1. Introduction 

The problem of a spin- 1/2 chain with anisotropic antiferromagnetic exchange interaction 
has been extensively studied [1] and constitutes one of best known examples of strongly 
correlated one-dimensional systems [2]. The XXZ model is integrable and exactly 
solvable by Bethe Ansatz [31 H] , which makes it possible to calculate exact ground state 
properties as well as thermodynamic quantities. At the same time, it exhibits a critical 
regime as a function of the anisotropy parameter, in which the system falls into the 
universality class of the Luttinger liquids. The long distance asymptotics of correlation 
functions can then be calculated by applying field theory methods. The combination 
of field theory and Bethe Ansatz has proved quite successful in explaining low energy 
properties of spin chain compounds such as Sr 2 Cu0 3 and KCuF 3 [5]. 

Recently, most of the interest in the XXZ model has turned to the study of 
dynamical correlation functions. The relevant quantities for spin chains are the 
dynamical structure factors S /J ' IJ '(q,u), \x = x,y,z, defined as the Fourier transform 
of the spin-spin correlation functions [6]. These are directly probed by inelastic neutron 
scattering experiments [TJ [8]. They are also probed indirectly by nuclear magnetic 
resonance [9] , since the spin lattice relaxation rate is proportional to the integral of the 
transverse structure factor over momentum [TDT [TT] . 

Even though one can use the Bethe Ansatz to construct the exact eigenstates, the 
evaluation of matrix elements, which still need to be summed up in order to obtain the 
correlation functions, turns out to be very complicated in general. In the last ten years 
significant progress has been made with the help of quantum group methods [12] . It is 
now possible to write down analytical expressions for the form factors for the class of two- 
spinon excitations for the Heisenberg chain (the isotropic point) at zero field [13 1 [HI [To] . 
as well as for four-spinon ones [T6[ [T71 [18] . No such expressions are available for general 
anisotropy in the gapless regime or for finite magnetic field, but in those cases the form 
factors can be expressed in terms of determinant formulas [T9], [201 IS] which can then 
be evaluated numerically for finite chains for two-particle states [22l l23l [24] or for the 
general multiparticle contributions throughout the Brillouin zone [25] |2~T]. 

From a field theory standpoint, dynamical correlations can be calculated fairly 
easily using bosonization [26]. However, this approach is only asymptotically exact in 
the limit of very low energies and relies on the approximation of linear dispersion for the 
elementary excitations. In some cases, the main features of a dynamical response depend 
on more detailed information about the excitation spectrum of the system at finite 
energies - namely the breaking of Lorenz invariance by band curvature effects. That 
poses a problem to the standard bosonization approach, in which nonlinear dispersion 
and interaction effects cannot be accommodated simultaneously. For that reason, a 
lot of effort has been put into understanding ID physics beyond the Luttinger model 

^ZII2HII2SII3DIEIII321I331I311I351I3S]. 

In particular, using the bosonization prescription one can relate the longitudinal 
dynamical structure factor S zz (q, u) at small momentum q to the spectral function of 
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the bosonic modes of the Luttinger model. In the linear dispersion approximation, the 
conventional answer is that S zz (q, oS) is a delta function peak at the energy carried by the 
noninteracting bosons [2|. As in the higher-dimensional counterparts, the broadening 
of the peak is a signature of a finite lifetime. The problem of calculating the actual 
lineshape of S zz (q, uj) at small q is thus related to the fundamental question of the 
decay of elementary excitations in ID. 

In the bosonization approach, interactions are included exactly, but band curvature 
effects must be treated perturbatively. All the difficulties stem from the fact that 
band curvature operators introduce interactions between the bosons and ruin the exact 
solvability of the Luttinger model. To make things worse, perturbation theory in those 
operators breaks down near the mass shell of the bosonic excitations j37j and no proper 
resummation scheme is known to date. The best alternative seems to be guided by 
the fermionic approach, which treats band curvature exactly but applies perturbation 
theory in the interaction [28] . 

In this paper we address this question using both bosonization and Bethe Ansatz. 
Our goal is to make predictions about S zz (q, u) that are nonperturbative in the 
interaction (i.e., anisotropy) parameter and are therefore valid in the entire gapless 
regime of the XXZ model (including the Heisenberg point). We focus on the finite field 
case, which in the bosonization approach is described by a simpler class of irrelevant 
operators. To go beyond the weakly interacting regime we can resort to the Bethe 
Ansatz equations in the thermodynamic limit to calculate the exact coupling constants 
of the low energy effective model. Our analysis is supported by another type of Bethe 
Ansatz based method, which calculates the exact form factors for finite chains. This 
provides a nontrivial consistency check of our results. 

The outline of the paper is as follows. In section we introduce the longitudinal 
dynamical structure factor for the XXZ model in a finite magnetic field and review the 
exact solution for the XX model. In section[3]we describe the effective bosonic model and 
explain how to fix the coupling constants of the irrelevant operators. Section H] provides 
a short description of the Bethe Ansatz framework which is relevant for our analysis. 
In section [5], we show how to obtain the broadening of S zz (q,u) in a finite magnetic 
field both from field theory and Bethe Ansatz and compare our formula with the exact 
form factors for finite chains. In section [6] we present a more detailed derivation of the 
high-frequency tail of S zz (q, u>) reported in [32]. The zero field case is briefly addressed 
in section [7J Finally, we check the sum rules and discuss the finite size scaling of the 
form factors in section [HJ 

2. XXZ model 

We consider the XXZ spin-1/2 chain in a magnetic field 



N 




(2.1) 
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Here, J is the exchange coupling, A is the anisotropy parameter, h is the magnetic 
field in units of J and N is the number of sites in the chain with periodic boundary 
conditions. We focus on the critical regime (given by — 1 < A < 1 for h = 0). We are 
interested in the longitudinal dynamical structure factor at zero temperature 

S zz (g,w) = Tf^2 e- lQU ~ f) / dte^ (S z (t) S], (0)> , (2.2) 

where q takes the discrete values q = 2irn/N, n e Z. It is instructive to write down the 
Lehmann representation for S zz (q,uj) 

S zz (q,u) = ^ Yl I <° I S l I a ) 1 2 6 ( w " E « + Egs ^ ' ( 2 ' 3 ) 

where S z = V. S z e~ tgj \a) is an eigenstate with energy i? Q and ^5 is the ground 
state energy. The matrix elements (0 \S Z \ a) are called form factors. We denote by 
F 2 = | <^ \ S* \ a) | 2 the transition probabilities that appear in (12.31) . For a finite system, 
S zz (q, uj) is a sum of delta function peaks at the energies of the eigenstates with fixed 
momentum q. In this sense, S zz (g, u) provides direct information about the excitation 
spectrum of the spin chain. In the thermodynamic limit N — »• oo, the spectrum is 
continuous and S zz (q, uj) becomes a smooth function of q and uj. Equation (12. 3p also 
implies that S zz (g, uj) is real and positive and can be expressed as a spectral function 

S* z (q,u) = -2Im X ret (q,cj), (2.4) 

for uj > 0. x ret (<?) u ) is the retarded spin-spin correlation function and can be obtained 
from the Matsubara correlation function 

i N r 13 

X(q,iu n ) = -jtJ2 / dre^(S z (r)S z ,(0)), (2.5) 

where (3 is the inverse temperature, by the analytical continuation iuj n — > uj + ie. 

It is well known that the one-dimensional XXZ model is equivalent to interacting 
spinless fermions on the lattice. The mapping is realized by the Wigner- Jordan 
transformation 

1 



3 3 2 



S+ -> {-1) J cje^, (2.6) 
Sj -> (-l) j Cj e- i7r ^, 

where c^- is the annihilation operator for fermions at site j, nj = c^c- and 0j = X^=i 
In terms of fermionic operators, the Hamiltonian (12.11) is written as 



TV 

- (ctc i+1 + /i.e. J -/ifctc.-- 



i\ / i 



+A[n i --jln i+1 2 



(2.7) 
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2.1. Exact solution for the XX model 

One case of special interest is the XX point A = 0, at which (12.71) reduces to a free 
fermion model [38] . As the free fermion point will serve as a guide for the resummation of 
the bosonic theory, we reproduce the solution in detail here. For A = the Hamiltonian 
(12. 7h can be easily diagonalized by introducing the operators in momentum space 

N 



Cp 



N . 



with p = 27m/ N, n e Z, for periodic boundary conditions. The free fermion Hamiltonian 
is then 

v 

where e p = — J (cos p + h) is the fermion dispersion. In the fermionic language, the 
dynamical structure factor reads 

i r+oo 

S zz {q,u>) = - J dte^(n q (t)n. q (0)) 
2tt 

= -jf I <° Kl «> I 2 6 (<" - E a + Eqs) , (2.10) 

a 

where n q = £V e'^nj = Y, p c\c p+q . 

We construct the ground state |0) by filling all the single-particle states up to the 
Fermi momentum kp. The latter is determined by the condition = 0, which gives 

kp = arccos (—h) — ir ( - + a ) , (2-H) 



2 

where a = (S"|) = (rij) — | is the magnetization per site. We can also describe the 
excited states in terms of particle-hole excitations created on the Fermi sea. The only 
nonvanishing form factors appearing in S zz (q, u) are those for excited states with only 
one particle-hole pair carrying total momentum q: \a) — c p+q c p |0). The form factors 
are simply 

(0 \S z q \a) = 6 (kp -\ p \)6(\p + q\- k F ) . (2.12) 

For a finite system there are qN/2n states with form factor 1, corresponding to different 
choices for the hole momentum p below the Fermi surface. In the limit iV — > oo, (12.101) 
reduces to the integral 

dpQ (k F - \p\) 9 (\p + q\- k F ) 5(lu- e p+q + e p) 

= 0(uj-uj l (g))9(uu (q)-u) ^ 
(dupjdp)^^ 

where u pq = e p+q — e p is the energy of the particle-hole pair and ujl (q) and ujjj (q) are the 
lower and upper thresholds of the two-particle spectrum, respectively. For the cosine 
dispersion, we have 

Lo pq = 2Jsin [p+ |^ sin|. (2.14) 
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Figure 1. Exact dynamical structure factor S zz (q, ui) for the free fermion point A = 0. 
For this graph we set a = —0.1 (kp = 2tt/5) and q = 7r/50. 



The expressions for the lower and upper thresholds depend on the proximity to half- 
filling (zero magnetic field). Here we shall restrict ourselves to finite field and small 
momentum \q\ -C kp. More precisely, we impose the condition 



\q\ < \2kp — 7r| = 2n\a 
For kp < 7r/2 (a < 0), we have 

(q) 



2 J sin — sin kp 

2 V 2 



,r. 9 • , \Q\ 
2J sm — sin \ kp H 

2 V 2 



(2.15) 

(2.16) 
(2.17) 



If fcp > 7r/2, the above expressions for ujp (q) and u>u (q) are exchanged. Hereafter we 
take k F < tt/2 and g > 0. It follows from ([2~T6D and fl2~T7l) that S zz for fixed g is 

finite within an energy interval of width 

5uj q = u>u(q) — u>L(q) = 4 J cos kp sin 2 J ^ (J cos kp)q 2 (2.18) 

for small g. In fact, we can calculate S* z2: (g, ou) explicitly using (12.131) . The result is 

9 (u - u L (g)) 9 (cuu (g) - lu) 



S zz (q,u) 



(2.19) 



(2 J sin |) -cu 2 

which is illustrated in figure [U Note that, although the form factors are constant, 
S zz (q,cu) is peaked at the upper threshold because of the larger density of states. The 
values of S zz (g, to) at the lower and upper thresholds are both finite 



S zz (g,w-> u L>u (g)) 



2 J sin | cos {kp =F ^ 



(2.20) 
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In the small-g limit, only excitations created around the Fermi surface contribute 
to S zz (q, uj). For this reason, a simplifying approach would be to expand the fermion 
dispersion around the Fermi points 

e R k L ^±v F k + ^ T ^ + ..., (2.21) 
K 2m 6 

where k = p =F kp for right (R) or left (L) movers, vp = J sin hp is the Fermi velocity, 

m = ( J cos kp)' 1 is the effective mass at the Fermi level and 7 = J sin kp. The free 

fermion Hamiltonian is then approximated by 

00 

#0 = Yl [ e * : cikRCkR ■ + e * : c ^ Cfci : ] ' (2 - 22) 

where CkR,p are the annhilation operators for fermions with momentum around ±k F , 
respectively, and : : denotes normal ordering with respect to the ground state. If we 
retain only the linear term in the expansion, Uk q turns out to be independent of k. This 
means that all particle-hole excitations are degenerate, and S zz (q, u>) is given by a single 
delta function peak at the corresponding energy lo = v F q 

S zz (q,u) = q5(uj -v F q) . (2.23) 

This is a direct consequence of the Lorentz invariance of the model with linear dispersion. 
In order to get the broadening of S zz (q, uj), we must account for the nonlinearity of the 
dispersion, i.e., band curvature at the Fermi level. If we keep the next (quadratic) term 
in e R,L , we find 

m f q 2 



S»{q,u,) = -e\^-\u,-v F q\). (2.24) 

We note that this flat distribution of spectral weight is a good approximation to the 
result in (12.181) and (12.191) in the limit q <C cot kp, in the sense that the difference 
between the values of S zz (q, uo) at the lower and upper thresholds is small compared 
to the average height of the peak (see figured]). This difference stems from the energy 
dependence of the density of states factor 1/ (dujp q / dp)]^^, which is recovered if we 
keep the fc 3 term in the dispersion. It is easy to verify that for q cot kp ('-/mq 1) 

AS ZZ = S zz (q, uu(q)) - S zz (q, u L (q)) « 7m 2 . (2.25) 

AS ZZ is q- independent, therefore AS ZZ /(m/q) ~ q vanishes as q — ► 0. This means that 
if we compare S zz (q, u) for different values of q - taking into account that Su q ~ q 2 and 
S zz {q,u) ~ 1/q inside the peak and rescaling the functions accordingly - the rescaled 
function becomes flatter as q — >• 0. On the other hand, the slope dS zz /du near the 
center of the peak diverges as q — > 0. 

The thresholds for the two-particle continuum, 

o 2 

wu,l (<?) « v F q ± j-, (2.26) 

are easy to interpret. For kp < n/2, the lower threshold corresponds to creating a hole 
at the state with momentum q below kp (a "deep hole") and placing the particle right 
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above the Fermi surface, whereas the upper one corresponds to the excitation composed 
of a "high-energy particle" at kp + q and a hole right at the Fermi surface [28] . 

Alternatively, we could have calculated the density- density correlation function, 
which for A = is given by the fermionic bubble 

dk 6(-k)9(k + q) 



2tt iuj 

Using the quadratic dispersion f 

m . I iuj 



- [UJ 



-uj 



X 



IUJ 



log 



Vpk + k 2 /2m, we find 
vpq + q 2 /2m\ 



- (UJ 



-uj 



(2.27) 



(2.28) 



2-rrq ° \iuj — vpq — q 2 /2m J 

The result ( 12.241) is then obtained by taking the imaginary part of x ret {<!■> u ) according 
to (El). 



3. Low energy effective Hamiltonian 

3.1. The free boson Hamiltonian 

For a general anisotropy A / 0, the Hamiltonian (12.71) describes interacting spinless 
fermions. The standard approach to study the low-energy (long-wavelength) limit of 
correlation functions of interacting one- dimensional systems is to use bosonization to 
map the problem to a free boson model - the Luttinger model pQ. This approach has the 
advantage of treating interactions exactly. As a first step, one introduces the fermionic 
field operators ipR,i (x) 



ip (x) 



D ikpx 



ipR (x) + e- ik ^ L (x) 



(3.1) 



+A 



*I>R,L 0) = —f= Y] c kR>L e 



(3.2) 



k=-A 



where L = N is the system size (we set the lattice spacing to 1) and A < n is a 
momentum cutoff. In the continuum limit, the kinetic energy part of the Hamiltonian 
in (I2.22p can be written as 



dx < : ip 



v F (-id x ) 



-ida.) 



2m 



R ■ 



v F (-id x ) + 



-id r 



2m 



(3.3) 



The 1/m term is usually dropped using the argument that it has a higher dimension 
and is irrelevant in the sense of the renormalization group. However, it introduces 
corrections to the Luttinger liquid fixed point which are associated with band curvature 
effects. Similarly, if we write the interaction term in (12.71) in the continuum limit, we 
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H, 



int 



AJ 



AJ 



dx : -0' (x) ijj (x) : : ijr (x + 1) -0 (x + 1) : 



cfe {p fi (x) p R (x + 1) + p L (x) p L (x + 1) 



+ p^ (x) p L (x + 1) + pl (x) p R (x + 1) 



+ 
+ 



«2fcF„/,t 



Vr ( X ) ( X ) WL ( X + X ) ^ii ( X + !) + ^- C - 
-i2fc F (2x+l)^ ^ ^ ^ + ^ ^ ( x + 1) + /i.e. 



(3.4) 

where pr,l = : yrl^rl '■• The ^ as ^ term corresponds to Umklapp scattering and is 
oscillating except at half-filling (where = 2ir). We will neglect that term for the 
finite field case, but will restore it in section [7] when we discuss the zero field case. 
We now use Abelian bosonization and write the fermion fields as 

1 



\Jlixa 



(3.5) 



where a ~ k F x is a short-distance cutoff and (j)R,L are the right and left components of 
a bosonic field d> and its dual field 9 



>R 







V2 ' 

^L + <j)R 



s/2 



(3.6) 
(3.7) 



which satisfy [0 (x) , d x '9 (x')] = i5 (x — x'). The density of right- and left-moving 
fermions can be shown to be related to the derivative of the bosonic fields 

Pr,l ~ =F— =d x (j) R , L , (3.8) 



so that 



n (x) 



1 



2tt 



cos I V47T0 — 2kpx 



(3.9) 



2 2vra 
Here we are interested in the uniform (small q) part of the fluctuation of ~ n(x), 
which is proportional to the derivative of the bosonic field <ft. Bosonizing the linear term 
in the kinetic energy (13. 3p . we find 

"L 



H, 



tin 





2 



dxiv F ( : ip R d x ip R : - : ^ L d x ^ L : 



dx [{d x (j) R f + (40l) 2 ] 



(3.10) 



The terms that appear in the interaction part are 

1 2 

Pr,l (x) Pr,l (x + 1) = — {d x (f) R>L ) , 

In 



PR (x) pL (X + 1) 



- fr-d x (j) R d x (j)L, 
Ztt 
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*Pr fa) i>L ( X ) i>\ ( X + 1) (X + 1) = 

cos(2fc F ) 2 . sin(2fc F ) 3 

- d^zj + - (d^ - <9 X .0 L ) + . . . , (3.11) 



2tt v ^ 3V2tt 

where we have set a = 1 (equal to the level spacing; see [2]). If we keep only the 
marginal operators (quadratic in d x (j)jn,), we get an exactly solvable model 



Hll = t / dx { f 1 + 2^;) [ ( ^ 0r)2 + ( ^ l)2] 



2 d x (j> L d x (S> R \, (3.12) 



7ryp 

where #2 = #4 = 2JA[1 — cos(2fcp)] = 4JAsin 2 /ci?. The Hamiltonian f)3.12p can be 
rewritten in the form 



hll= ij dx vK ( dx °y + ^ ( 9x 



(3.13) 



where w (the renormalized velocity) and K (the Luttinger parameter) are given by 

2A 




7T 



sinfcp , (3.14) 



K = < / i , ^ « 1 - V 8infc F . (3.15) 

2rrvp 2ttvf 

Expressions (I3.14p and (13.151) are approximations valid in the limit A < 1. The 
Luttinger model describes free bosons that propagate with velocity v and is the correct 
low energy fixed point for the XXZ chain for any value of A and h in the gap less regime. 
However, the correct values of v and K for finite A must be obtained by comparison 
with the exact Bethe Ansatz (BA) solution. In the case h = 0, the BA equations can 
be solved analytically and yield 



v(A,h = 0) =^^^, (3.16) 
2 arccos A 

K (A, h = 0) = — — . (3.17) 

v ' ; 2 (vr - arccos A) v ; 

There are also analytical expressions for h « and h close to the critical field [39]. For 
arbitrary fields, one has to solve the BA equations numerically in order to get the exact 
v and K. 

The Luttinger parameter in the Hamiltonian (13.131) can be absorbed by performing a 
canonical transformation that rescales the fields in the form — > \^Kcj) and 9 — > 6/\/K~. 
Hll then reads 

R LL = \j dx [{d x 6f + (d x <j>) 2 ] . (3.18) 
We can also define the right and left components of these rescaled bosonic fields by 
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These are related to (f)R t L by a Bogoliubov transformation. An explicit mode expansion 
(neglecting zero mode operators) is 

(x, r) = 7^ [af L e-^i*) + ^t^-)] ? (3 20) 

<2>0 

where a^ ,L are bosonic operators obeying [a^' L , a^,' L ^} = 5 qq > and q = 2nn/L, n = 
1,2, . . ., for periodic boundary conditions. The Hamiltonian (13.181) is then diagonal in 
the boson operators 

Hll = [afaf + a^a L q } . (3.21) 

g>0 

We can calculate the propagators for the free fields d x (p^L from the mode expansion 
in (13.201) . In real space, for L —> oo and zero temperature {(3 — > oo), the propagators 
read 

D ( Sl t) = (d xVR ,L (x, t) d x p R ,L (0, 0)) = ~ (3.22) 

Ztc [vt =f IX) 

In momentum space, 

D® L {q,iu n ) = - f dxe-^ I ' dr e^D^ L (x,r) 
Jo Jo 

= ±q . (3.23) 

iu n =F vq 

In order to calculate the dynamical structure factor defined in ( 12. 2ft , we express the 
fluctuation of the spin operator in terms of the bosonic field <fi. From (12.61) and (13. 9p . 
we have 

SI ~ \f^dA. (3.24) 



In the continuum limit, 



l r /3 



X(q,iw n ) = -— / tee-*** I dr e^ T {d x tf> (x, r) d x ^> (0, 0)) 

= ^D<®(q,i(j n ), (3.25) 
where D^(q,iu n ) is the free boson propagator (for the d x (f> field) 

L><°) (g, iu) = 4 0) (g, iu,) + Z>f (g, iu;) = - 2 W 2 . (3.26) 

(iw) - (vg) 

It follows that the retarded correlation function is 

1 1 



ret l \ K 1 



(3.27) 

2tx \_uj — vq + irj lo + vq + ir\ 
Finally, using (12 .4p . the dynamical structure factor for the free boson model is (g > 0) 

S zz (q,u) = Kq5{u-vq). (3.28) 

The result in (13.281) is analogous to (I2.23p . Since the Luttinger model exhibits Lorentz 
invariance, S zz (g, u) is a delta function peak at the energy carried by the single 
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boson with momentum q. This solution should be asymptotically exact in the limit 
q — > 0, which means that any corrections to it must be suppressed by higher powers 
of momentum. However, the free boson result misses many of the features that the 
complete solution must have. For example, the exact solution for the XX point suggests 
a broadening of the delta peak with a width Su q ~ q 2 . Like in that case, it is necessary 
to incorporate information about band curvature at the Fermi level by keeping the 
quadratic term in the fermion dispersion in order to get a finite width for S zz (q,uj). 
As we shall discuss in the next section, the problem is that such a term is mapped 
via bosonization onto a boson-boson interaction term. Even though the interaction 
term is irrelevant, finite-order perturbation theory in these operators leads to a singular 
frequency dependence close to uj = vq. It turns out that broadening the delta function 
peak within a field theory approach is a not an easy task. A complete solution that 
recovers the scaling 5uj q ~ q 2 requires summing an infinite series of diagrams, as we 
will point out in section [51 Another feature expected for S zz (q, uj) when A ^ is a 
high-frequency tail associated with multiple particle-hole excitations. This tail can be 
calculated in the region 5uj q <C oj — v q <C J by lowest-order perturbation theory in 
the fermionic interaction (oc A) starting from a model of free fermions with quadratic 
dispersion [27] . In section [6] we obtain this result by including fermionic interactions 
exactly (finite A) and doing perturbation theory in the band curvature terms. 

3.2. Irrelevant operators 

In order to go beyond the Luttinger model, we need to treat the irrelevant operators 
that break Lorenz invariance. There are two sources of such terms: band curvature 
terms, which are quadratic in fermions but involve higher derivatives, and irrelevant 
interaction terms [33]. The first type appeared in (13.31) and corresponds to the k 2 term 
in the expansion of the fermion dispersion 



We derive the bosonized version of a general band curvature term in the following way 
(see [10]). We define the operator 




(3.29) 





(3.30) 



where ... is a total derivative. Organizing by powers of e, we can write 




(3.31) 



n=0 



where 




(3.32) 
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(3.33) 



\j2na \[L 

where 0^ are the creation and annihilation parts of <pR (x) = 0^ (x) + 0^ (x) and we 
have used the identity e A+B = e A e B e~^ A,B ^ 2 with 



R \ X ) > 'rR 



(y)] 



2tt 



log 



2m 



[X — y + iac) 



(3.34) 



for large L. Then we express F (x, e) in terms of the bosonic fields 

F f Xj A = £ e i^^(x+e/2) e i^0-(a!+e/2) e -i^0+(x-e/2) e -i^0^(x-e/2)^3_3 5 ^ 

L 

After normal ordering the operators, we can do the expansion in e (dropping the normal 
ordering sign) 

tJ>r ( x + |) $r ( x ~ \ 



exp { i\j2Tx 

2nie P 



\<f>R (a 



2vrie£! ^ FT rnA V2 



n 

3=1,3,- 



dj(pR 



(3.36) 



From (I3.3ip and the coefficient of the e n term in (I3.36p . we have 



x 



2 n+1 2m 



E 



(2-s/27r«) 



n 



(3.37) 



where the m/s obey the constraint Y2jj m j = n + 1. In particular, for n = 2 the sum 
in (I3.37P contains only two terms (either mi = 3, ms = or mi = 0, 7713 = 1). We get 



f {2) (x) = vJAV R 



'2tt ,„ , , 3 



1 



9|0R. 



(3.38) 



3 v 12V27T 

The last term is a total derivative and can be omitted from the Hamiltonian. Similar 
expressions for the left-moving field 0l are obtained straightforwardly by using the 
symmetry under the parity transformation x — > —x, R — > L. The bosonized version of 
the band curvature terms in (I3.29P is then 

^7T 



5H 



be 



6m 



(3.39) 



We now rewrite 57ib c in terms of the right and left components of the rescaled field 
). Using (EBD and <K7h . 



sn 



be 



6m 



d j + d x <, 
V2 



~\ 3 



d x 9 - d x 
V2 



I dx [3 (d x 6) 2 d x( j> + K 2 («9 X .0) 3 ] . 



(3.40) 
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Finally, using (13.191) . we get (in accordance with [35J) 

SH bc = ^- [(d x <p L ) - (d x <p R ) \ 



[(Ob^l) 0^ - (S^) d^zj • (3.41) 



6 4m 

Besides 57ib c , we need to include the irrelevant operators which arise from the 
expansion of the fermionic interaction in the lattice spacing, as we encountered in ( 13. lip . 
In terms of <p Rt L, this contribution reads 

JAK 3 / 2 

5H int = - sin(2fc F ) { [(«9^ L ) 3 - (d x <p R ) 3 ] 
3v 27T 

-3 [(d^) 2 fie^a - (d x ip R ) 2 d x ip L ] } . (3.42) 

Combining (I3.4ip and (I3.42p . we can write the irrelevant operators in the most general 
form 

5H = ^ j dx {m_ [(d xVL ) 3 - (d x cp R f] 

+?7+ [{d x (p L ) 2 d x <f R - (d x <f R ) 2 d x ip L ] } . (3.43) 
To first order in A, the coupling constants r)± are given by 

w — f 1 + — sinA; F J , (3.44) 



m \ 7r 



3A 

77+ ~ sinfcp. (3.45) 

7rm 

The perturbation 5H in (I3.43P might as well have been introduced phenomenolog- 
ically in the effective Hamiltonian. In fact, the dimension-three operators (d x (p R! i,) 3 are 
the leading irrelevant operators that are allowed by symmetry. They obey the parity 
symmetry cfL — > f R , x —x, but not spin reversal (or particle-hole) f R ,L —<Pr,l, 
which is absent for ft ^ 0. Such terms give rise to three-legged interaction vertices 
which scale with powers of the momenta of the scattered bosons (figure [2]). They are 
responsible, for example, for corrections to the long distance asymptotics of the cor- 
relation functions [10]. Note that as A — ► (K — > 1), m_ — > 1/m while r] + vanishes 
because there is no mixing between right and left movers at the free fermion point. 
Moreover, the weak coupling expressions predict that both r]_ and r] + vanish in the 
limit h — >• (m , oo), in which particle- hole symmetry is recovered. (See, however, 
figure [15] below.) For h = the leading irrelevant operators are the dimension- four 
operators (d x ip R> L) , {d x ip R ) 2 (c^y?^) 2 and the umklapp interaction cos(4\A"-ft"0), which 
becomes nonoscillating [32] . 

The condition that a general model of the form Hll + SH be unitarily equivalent 
to free fermions up to dimension-four operators [33] amounts to imposing that the 
Bogoliubov transformation that diagonalizes Hll in the R/L basis also diagonalizes the 
cubic operators in 5H. In our notation, this condition is expressed as r] + = 0. That 
condition is not satisfied by the XXZ model except for the trivial case A = 0. However, 
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Figure 2. Interaction vertices in the low energy effective Hamiltonian. The solid 
(dashed) lines represent propagators for right- (left-) moving bosons D ( °\D { ^). 



the contributions from this extra (i.e., not present for free fermions) dimension-three 
operator to S zz (q,u) are of 0(t]^), as we will discuss in section El 

Similarly to what happens for v and K, (13.441) and (13.451) should be regarded as 
weak-coupling expressions. Again we can use the fact that the XXZ model is integrable 
and obtain the exact (renormalized) values of t]± by comparison with Bethe Ansatz. 
In section 13.31 we will discuss how to fix these coupling constants in order to obtain a 
parameter-free theory. 



3.3. Determination of the renormalized coupling constants 



As mentioned in section 13. 2\ the renormalized parameters r]± can be determined by 
comparison with exact Bethe Ansatz results for infinite length. We will proceed by 
analogy with the calculation for the zero- field case in [H]. One difficulty is that there 
are no analytical solutions of the Bethe Ansatz equations for finite fields, so we must be 
satisfied with a numerical evaluation of the parameters. In the following, we will relate 
rj± to the coefficients of the expansion of v and K as functions of the magnetic field, by 
comparing the corrections to the free boson result for the free energy calculated in two 
different ways. 

Let us consider the response to a small variation in the magnetic field around a 
finite value h . In the limit 5h = h — li < 1, such response is well described by the 
Luttinger model 



H = dxl- \{d x 9) 2 + (c\0) 2 l - J6hJ-d x <f> } , (3.46) 




where v (h) and K (h) are known exactly from the Bethe Ansatz equations. For ho = 0, 
the cutoff-independent terms of the free energy density according to field theory read 

/( fc = )„-£- *-(.H k )>, (3.47) 

where v and K are given by (13 . 1 6[) and (I3.17p . respectively. The magnetic susceptibility 
at zero temperature is x = —J 2 {d 2 f /dh 2 )\ T=0 = K/ttv, which is the familiar free 
boson result. For finite field ho =/ 0, the free energy assumes some general form 

/(^°)~-^ ~C(h), (3.48) 
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and the T = susceptibility is obtained by 



X = — « 



1 (d 2 f 



J 2 \dh 2 



1_ f d 2 C 
■P \d(Shf 



h,T=0 " W V Jlb J / h,T=0 

where the last identity holds for any Luttinger liquid. 

We would like to calculate the corrections to / and x that involve higher powers 
of the perturbation Sh. Our first approach is to assume that the field dependence is 
already completely contained in the definitions oiv(h) and K(h), so that we can employ 
the expansion 

v (h) = v (ho) [1 + aSh + (5h 2 )] , (3.50) 
K(h) = K(ho) [l + bSh + O (5h 2 )] , (3.51) 

where the coefficients a and b can be extracted from the exact v and K by linearizing 
the field dependence around h = h$. Consequently, the lowest-order correction to the 
free boson susceptibility around h = ho is 

X = ^Q[l-(°-b)6h + 0(6h 2 )]. (3.52) 
irv[ho) 

Likewise, the free energy at finite temperature must contain a term of the form 
ii5hT 2 

Sf ~ a — — — , (3.53) 
6v(h ) 

due to the field dependence of the velocity. Both a and b depend on ho and the anisotropy 
A. As an example, at the XX point, K = 1 for any value of the field, therefore 
6(A = 0, ho) = 0. From fl2TTTD . we have 

v F = J sin k F = jVT^h 2 w v F (h ) Sh + O (5h 2 ) , (3.54) 

v F (ho) 

so that we get 

^ a n t, \ j2h ° coskp a ^\ 

a A = 0,/i = 27TT = • 2 , • 3 - 55 

v F {ho) sm k F 

In our second approach, we take v = v(h ) and K = K(ho) to be fixed and assume 
that the corrections to the free boson result are generated by the irrelevant operators. 
We consider the effective Hamiltonian H = Hll + SH, with SH defined in (13.431) . An 
equivalent Lagrangian formulation in imaginary time is 

C = C + SC, (3.56) 

Co = + V - (<9,0) 2 - JShMd x <P, (3.57) 

2v 2 V vr 

6C = ~ i^r (<^) 2 d ^ + ^ (<^) 3 + W±) . ( 3 - 58 ) 
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where A = 3r/_ +r/ + and B = rj- —rj + . We shift the field by 4> — > <p+ ^py^s to absorb 
the term linear in d x <p and get 

£o .(M! + W + «, (3 .59) 
2v 2 2irv 

BK 3 ' 2 (J5h) 3 , , r , 

H ^ h odd powers or 0. (3.60) 

67TU d 

We then calculate the free energy density from the partition function 

Z = J V(f> exp |- J dr dx(C + 5C)^, (3.61) 

rp rp p/3 



f = __ i nZ ~ fo + _ rfr rfx 9 (3.62) 

where fo reproduces the free boson result 

and (5C) is the expectation value of 5C calculated with the unperturbed Hamiltonian. 
In order to compute (SC), we need the finite temperature propagators 

(d x <p(x + e) d x (j>{x)) = - \ (d T <f) {x + e) d T <p {x)) 

1 (nT/vf ^ 
2tx sinh 2 (irTe/v) 

Now we use the expansion sinh -2 (jrTe/v) ~ (v/irTe) 2 — 1/3 for e — ► and drop the 
cutoff-dependent terms in 5f. The reason is that the latter simply renormalize the 
corresponding terms in fo and have already been accounted for in the renormalization 
of v and K. The correction to the free energy to first order in A and B becomes 



rp i>P i-L 

5f = - dr dx (5C) 
L Jo Jo 



it \f~KJbh T 2 K 3 ' 2 (J8hf 



{A + W) 36, +B ^^- (3 ' 65 » 



The susceptibility obtained from fo + 5f is 



1 d 2 (f + 5f) 

X T9 rwr,x2 



JL_ B ™ +0( ^). (3 . 66) 



1XV 1XV 



J 2 d(5hy 

Comparing with the expression (I3.52p . we can identify 

a-b = — —B = — — (??_ -rj+). (3.67) 

Besides, from the ShT 2 term in (I3.53P and (13.651) . we have 

VKJ , t . Vkj , . 

a = ^r( A + w ) = ^r - n+) ■ ( 3 - 68 ) 
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Finally, combining (I3.67P and (I3.68P and writing a = v 1 dv/dh and b = K 1 dK/dh, 
we find the formulas first presented in [32] 

v dv v 2 dK 
3v 2 dK 

Jv+ = 2K^m- (3 - 70) 

The above relations allow us to calculate the renormalized values of r)± once we have 
the field dependence of v and K. Notice that r] + oc dK/dh and as expected vanishes 
at the XX point. On the other hand, r\- remains finite at A = because dv/dh ^ 
and we recover rj_ = (vp / J)dvp/dh = J cos hp = m -1 . It is also possible to check the 
validity of (13. 69ft and ( 13.701) explicitly in the weak coupling limit, using the expressions 
for t)(A<l,/i) and K (A <C 1, h) in ( 13.141) and ( 13. 15ft as well as the weak coupling 
expressions for r]± in (13.441) and (13.451) . 



4. Bethe Ansatz solution 



Although the Bethe Ansatz is first and foremost a method for calculating the energy 
levels of an exactly solvable model (readers who are unfamiliar with the subject 
are invited to consult standard textbooks, for example [12J [391 H3])> recent progress 
stemming from the Algebraic Bethe Ansatz means that we can now use it to make 
many nontrivial statements about dynamical quantities. Assuming that certain specific 
families of excited states carry the dominant part of the structure factor, we can delimit 
the energy and momentum continua where we expect most of the correlation weight 
to be found, and provide the specific lineshape of the structure factor both within this 
interval, and further up within the higher-energy tail. We start here by introducing 
the important aspects of the Bethe Ansatz which we will make use of later on when 
studying the correspondence with field theory results. 



4-1. Bethe Ansatz setup and fundamental equations 

As is well-known, an eigenbasis for the XXZ chain (12.11) on N sites is obtained from 
the Bethe Ansatz PS], 

^a/(Ji, ...,3m) = ^(-l)l p 'e'^i^ 3 »-^i< <KM#p a M). (4.1) 
p 

Here, M < N/2 represents the number of overturned spins, starting from the reference 
state |0) = ®f = i\ t)i (i- e - the state with all spins pointing upwards in the z direction). 
The total magnetization of the system along the z axis, S^ ot = Na = y — M is conserved 
by the Hamiltonian. P represents a permutation of the integers {1, ...,M} and ji are 
the lattice coordinates. The quasi- momenta k are parametrized in terms of rapidities A, 

e,t = Si ta + r/2V A = cos <- < 4 - 2 > 

sinh (A — 



Dynamical structure factor at small q for the XXZ spin- 1/2 chain 



19 



such that the two-particle scattering phase shift becomes a function of the rapidity 
difference only, 4>(k a , = <fti(\ a — A&) with <pi defined below. An individual eigenstate is 
thus fully characterized by a set of rapidities {A}, satisfying the quantization conditions 
(Bethe equations) obtained by requiring periodicity of the Bethe wavefunction (14.11) : 

1 M I 
- ]v E M A i - A *) = 27r ]f 3 = 1,-, M, (4.3) 

k=l 

in which Ij are half-odd integers for N — M even and integers for N — M odd, and where 
we have defined the functions 

n (A) = 2 arctan f . (4.4) 



_tan«/2) 

The energy and momentum of an eigenstate are simple functions of its rapidities, 
E = — tt J sin ( ai(\j) — hS^ 



J tot, 



P = nM _ 0i (A,-) = ttM - — J] (4.5) 



in which 



1 d , 1 sin(nC) 

27T gsA 7r cosh(2A) — cos(n^) 

Each solution of the set of coupled nonlinear equations (14. 3 j) for sets of non- 
coincident rapidities represents an eigenstate (if two rapidities coincide, the Bethe 
wavefunction (14. ip formally vanishes). The space of solutions is not restricted to 
real rapidities: it has been known since Bethe's original paper that there exist 
solutions having complex rapidities ('string' states), representing bound states of 
magnons. In fact, obtaining all wavefunctions from solutions to the Bethe equations (or 
degenerations thereof) remains to this day an open problem in the theory of integrable 
models. It is however possible to construct the vast majority of eigenstates using 
this procedure, allowing to obtain reliable results for thermodynamic quantities and 
correlation functions. In all our considerations in the present paper, we can and will 
restrict ourselves to real solutions to the Bethe equations. 



4-2. Ground state and excitations 

The simplest state to construct is the ground state, which is obtained by setting the 
quantum numbers Ij to (we consider iV even from now on for simplicity) 

/f = -^+J, J = 1,-,M. (4.7) 

The simplest excited states which can be constructed at finite magnetic field are 
obtained by introducing particle-hole excitations on the ground-state quantum number 
distributions, see figure ([3]). Since we limit ourselves to real solutions to the Bethe 
equations, we require \Xj\ < oo and thus |/| < 1^, where, from (14.31) . 

N-M N (a o\ 
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Figure 3. Representation of various quantum number configurations: a black (empty) 
circle represents an occupied (unoccupied) allowable quantum number (which here are 
taken to be half-odd integers; the figure is centered on zero). The top set represents 
the ground state configuration, whereas the second and third from top represent two- 
particle excitations at different momenta, obtained by creating a particle-hole pair on 
the ground-state configuration. The bottom set is for a four-particle state obtained 
from two particle-hole pairs. The dotted line delimits the ground state interval, whereas 
the solid lines delimit the quantum numbers for which real solutions to the Bethc 
equations can be obtained in this illustrative case (see main text). 

The momentum of an excited state is simply given by the left-displacement of 
the quantum numbers with respect to those in the ground state, q = jffyl, where 
5il = X]j(-^ G5 ~~ h)- At a given fixed (small) momentum, we can thus construct qN/27r 
two-particle states by shifting the particle and hole quantum numbers, leaving their 
difference fixed. Since the energy of these two-particle states at fixed momentum are non- 
degenerate, this defines a two-particle continuum whose characteristics will be studied 
later. Higher-particle states can be similarly constructed and counted. 

The restriction to real rapidities and a single particle-hole pair therefore means 
that our subsequent arguments will apply only to the region q < Min(2fcp, koo) where 
kp = ii^r and k^ is given by the maximal displacement of the outermost quantum 
number, k^ = 27r / °°~ A f / ^ 2 . We can thus write our restriction as 

q< Min{7r(l-2o-),2o-(7r-C)} (4.9) 

in terms of the magnetization, noting in particular that the window of validity of our 
arguments vanishes in the case of zero magnetic field. 

For a finite chain with N sites and M overturned spins, the Hilbert space is finite, 
and therefore so is the sum over intermediate states in the Lehmann representation for 
the structure factor (12.31) . Each intermediate state is obtained by solving the Bethe 
equations, the space of states being reconstructed by spanning through the sets of 
allowable quantum numbers. The form factor of a local spin operator between the 
ground state and a particular excited state is obtained from the Algebraic Bethe Ansatz 
as a determinant of a matrix depending only on the rapidities of the eigenstates involved 
[T9l [20] even in the case of string states with complex rapidities [2T]. This enables to 
obtain extremely accurate results on the full dynamical spin-spin correlation functions 
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in integrable Heisenberg chains [23121]. We will make use of this method in what follows 
to compare results from the Bethe Ansatz to field theory predictions for the structure 
factor at small momentum. 

5. Width of the on-shell peak 

Linearizing the dispersion around the Fermi points is a key step for the bosonization 
technique. By doing so all the particle-hole excitations with same momentum q <C hp 
become exactly degenerate and one can associate a particular linear combination with a 
single-boson state [3D]. In this approximation, the single boson state \b) = |0) is the 
only state that couples to the ground state via S q . The associated weight in S zz (q, uj) 



However, as we will see in section l5~3l the exact eigenstates in the Bethe Ansatz solution, 
whose energies are given by (14.51) . are nondegenerate. In fact, most of the above spectral 
weight is shared by qN/27r two-particle states whose energies are spread around uj = vq. 
This is reminiscent of the exact solution for the free fermion point in section 12.11 In 
the bosonic picture, on the other hand, the broadening 5uj q is related to a finite lifetime 
for the bosons of the Luttinger model. Once band curvature is introduced via the 
irrelevant operators in (13.431) . the single boson is allowed to decay and the coupling 
to the multiboson states lifts the previous degeneracy. The fact that the irrelevant 
operators have the same scaling dimension as in the noninteracting case suggests that 
for A ^ the width should also vanish as q 2 in the limit q — > 0. In this section we 
argue in favor of a q 2 scaling for 8uj q for all values of A in the gapless regime, as long 
as r\- 0, based on two different approaches. First, we explain how the expansion 
of the bosonic diagrams in the interaction vertex r)-, neglecting rj + , coincides with the 
expansion of the free fermion result d 2 . 2 8 [) in powers of 1/m. rj_ is then interpreted 
as a renormalized inverse mass, in the sense that the width of the peak for A ^ is 
given by 5uj q = \rj_\q 2 . Second, we derive from the Bethe Ansatz equations an analytical 
expression for the width of the two-particle continuum at finite fields and show that it 
coincides with the field theory prediction for the width of S zz (q, uj). Finally, we confirm 
these results directly by analyzing the numerical form factors calculated for finite chains 
of lengths up to 7000 sites. 

5.1. Width from field theory 

We saw that the width 5uj q is well defined for the free fermion point, in which case 
S zz (q, uj) has sharp lower and upper thresholds uj LtU (q). For the interacting case, 
S zz (q, uj) still vanishes below some finite lower threshold ujl{q) at zero temperature 
due to simple kinematic constraints. However, the on-shell peak has to match a high- 
frequency tail somewhere around uj\j (q), hence the meaning of an upper threshold is no 
longer clear. 



is given by 




(5.1) 
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In their solution for weakly interacting spinless fermions, Pustilnik et al. [2S] found 
that ujjj has to be interpreted as the energy at which the peak joins the high frequency 
tail by approaching a finite value with an infinite slope. Although it is actually possible 
that the singularity at uiu (q) get smoothed out if one treats the decay of the "high- 
energy electron" for a general model [36] , the singularity may be protected in integrable 
models such as the XXZ model. 

Of course the situation is a lot simpler for models with no high-frequency tail, where 
the dynamical structure factor is finite only within the interval 0Ji{q) < uo < uojj(q). 
Such is the case for the Calogero-Sutherland model [H]. The absence of a tail for S (q, uj) 
in the Calogero-Sutherland model can be attributed to the remarkable property that 
the quasiparticles are all right movers [29] . As we will discuss in section |6l the t] + term 
that mixes R and L in our low energy effective Hamiltonian (figure [2]) is responsible 
for the high-frequency tail for h ^ because it allows for intermediate states with two 
bosons moving in opposite directions, thus carrying small momentum and high energy 
uj ^> vq. 

In contrast, the rj_ interaction has matrix elements between multiboson states which 
contain only right movers. All these states have uj ~ vq. Therefore 7/_ must be related 
to the broadening of the on-shell peak. It has already been pointed out in [33] that 
the model with 7] + = is equivalent to free fermions up to irrelevant operators with 
dimension four and higher. For this case one can write down an approximate expression 
for the dynamical structure factor which misses more subtle features in the lineshape 
(e.g., the power law singularities at the thresholds) but accounts for the renormalization 
of the width due to interactions. Even for models with nonzero such as the XXZ 
model in the entire gapless regime, it is reasonable to expect that 5uj q , if well defined, 
will depend primarily on the interaction between excitations created around the same 
Fermi point. For that reason, we will neglect the r] + interaction in an attempt to derive 
an expression for the width of S zz (q, uj) from the bosonic Hamiltonian. In the following 
we apply perturbation theory in rj_ up to fourth order and show that it recovers the 
expansion of the logarithm for the density-density correlation function. This fact has 
already been noticed in [34"[, [35] up to 0(r/^). However, irrelevant interaction terms 
such as (13. 42ft were neglected in [34], [35] . Such terms are crucial to obtain the correct 
effective inverse mass, since the correction of first order in the fermionic interaction A 
stems from (13.421) . 

For rj + = 0, the Hamiltoninan Hn + SH decouples into right and left movers. For 
excitations with q > 0, we can consider only right movers and work with 

V 2 V2"?T 3 

Hr = j ( 9 x<Pr) ^— V- (9x<fiR) ■ (5.2) 

The first attempt to broaden the delta function peak in S zz (q, uj) would be to calculate 
the corrections to the propagator 

TV- j-+oo rj3 

x ( q ,iuj) = -— / dxe-w / dre^ T (T T d x ip R {x,r) 0^(0,0)), (5.3) 

^ J-oo Jo 
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by using perturbation theory in the cubic term. Unfortunately, any finite order 
perturbation theory in r)_ breaks down near uj « vq. Even using the Born 
approximation, which sums an infinite series but not all the diagrams, one finds that 
the self-energy to 0(rfi) is divergent: Im £ (q, uj) ~ 5 (uj — vq) [37]. This is actually not 
surprising if we look at the exact solution for the free fermion point. Expanding the 
positive- frequency part of (12.281) in powers of 1/m, we get 

2 i / 9 / \ 4 



q 



2ttw 



2w 



q 2 /m 



2w 



(5.4) 



where w = iuj — vpq. Stricly speaking, such expansion is valid only for uj — vpq ^> q 2 /2m. 
For uj ~ Vpq, the expansion in band curvature produces increasingly singular terms that 
need to be summed up to produce the finite result in (12.281) . 

In any case, it is legitimate to examine the expansion of bosonic diagrams and ask 
whether it can at least reproduce the free fermion result. We use the bare propagator 



D 



(0) 
R 



x, t) = (T T d x tp R (x, r) d x (fR (0, 0)) in (13.221) . with Fourier transform 
Df (q, kj) 



(5.5) 



to calculate the expansion of xr (<L ^) U P to O (t/I), as represented in figure HI The 
zeroth-order result is simply the same as in (13.271) 

Kq 



The 0(r}i 



X {0) (q^)- 
correction is 



2ttw 



(5.6) 



X 



(2) 



lUJ) 



K 
2^ 



d xe 



-iqx+iuir 



d 2 Xi / d 2 x 2 



2tx 



x 



(T T d x <p R (x) [d x ip R (l)] 3 [d x ip R (2)] 3 d x ip R (0)) 



A' 

2^ L 



D 



(0) 
R 



lUJ] 



U RR (q, iuj) , 



(5.7) 



where H RR (q, iw) is the bubble with two right-moving bosons 

U RR (q,iu) = -tttiI Cf- ( + ™ d ^Df{kMDf{q-k 



ri 2 _q z 
T2^' 



{5.1 



Note that H RR is singular at uj = vq, which prevents us from treating it as a self-energy. 
The origin of the singularity is that the two right-moving bosons in the intermediate 
state always carry energy uj = vq, no matter how the momentum is distributed between 
the pair. Substituting (15.81) back into (15.71) . we get 

<2) («,M=#^(^y. (5.9) 



r 



27TW 12 



w 
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X = 




Figure 4. Perturbative diagrams up to fourth order in r\- 



To O (jit), there are three topologically distinct diagrams (figure HJ), which give the 
following contributions 

,2' 



ifg 1 



27TM7 144 

ifg 1 
2?™ 504 
ifg 1 
2mv 280 



^-3 



(5.10) 



The coefficients for each diagram are nontrivial and result from both combinatorial 
factors and integration over internal momenta (recall that the interaction vertex is 
momentum-dependent because of the derivatives in (13.43}) ) . Remarkably, all the fourth- 
order diagrams have the same q and to dependence with comparable amplitudes. We 
are not allowed to drop any of them and there is no justification for the use of a self- 



consistent Born approximation, for example 
up with the expansion 



Putting all the terms together, we end 



Kq 
2ttw 



?7_g 
~2w 



1 

+ 5 



rj-q 



+ 



(5.11) 



which is analogous to (15.41) with the replacements 1/m — *■ r\- , Vf — * v and an extra factor 
of K. This proves that the expansion of bosonic diagrams reproduces the expansion of 
the free fermion result up to fourth order in 1/m. Since there is no simple way to predict 
the prefactors of each diagram, all we can do is to check this correspondence order by 
order in perturbation theory. However, if we believe that the bosonic theory reproduces 
the free fermion result to all orders in rj_, we must conclude that in the interacting case 
the series in (15.111) sums up to give the result 



X 



K 



IUJ] 



27cr)_q 



log 



iuj — vq + ?7_g 2 /2 
iuj — vq — r]_q 2 /2 



from which we obtain 

S zz (q,u) 



K 

\v-\q 



-e 



r)-\q 



\uj — vq 



This result predicts that S zz (q, to) is finite and flat within an interval of width 



5u n 



\V-\Q 



(5.12) 



(5.13) 



(5.14) 
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Figure 5. Lincshape in the approximation with the r;_ interaction only (solid line). 
The dotted line illustrates the expected true lineshape for small A (see section [5731 ). 



This lineshape (illustrated in figure [5]) is the exact one for the case of free fermions with 
quadratic dispersion. The reason is simple: because the bosonization of the operator 
~ k 2 c\ck only generates the r\_ term, one could invert the problem and refermionize the 
Hamiltonian (15.21) to an effective free fermion model with inverse mass r\-. In a more 
general model, more irrelevant operator have to be added to the effective Hamiltonian to 
reproduce details of the lineshape that are higher order in q. For example, we expect the 
power-law singularities present at ujl,u f° r A ^ [28] to be associated with dimension- 
four operators such as (d 2 ipn) 2 and (c^y?^) 4 (with corrections of 0{rj 2 [ _), see section [7]). 
This means that if we write 



Q 



S zz (q iU ) = ^-f[q 



lu — vq 



(5.15) 

J i v UUJ 1 

the rescaled function f(q,x) approaches the flat distribution of figure [5] in the limit 
q — > 0. Finally, we note that this approximate solution yields the same sum rules as the 
free boson result 



duS zz (q,u) =Kq, 

duujS zz (q,u) = vKq 2 , 
and also the magnetic susceptibility 

X = x(g = 0) = lim- r^S zz (q,u) 

q^O TV Jq LU 

independent of the value of ??_. 



K 



TTV 



(5.16) 
(5.17) 



(5.18) 



5.2. Width from Bethe Ansatz 

The purpose of this section is to provide an analytical derivation of the quadratic width 
formula f)5.14p . making use of standard methods associated to the thermodynamic Bethe 
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Ansatz, and assuming that single particle-hole type excitations in the Bethe eigenstates 
basis carry the most important part of the structure factor. We first set our notations 
and underline certain characteristics of the ground state of the infinite chain in a field 
which will prove to be useful for our purposes. We then discuss particle-hole excitations 
in the thermodynamic limit, and obtain a relationship giving the width in terms of 
solutions of integral equations, simplifying to the conjectured field theory result in the 
small momentum limit. 

Let us begin by taking the thermodynamic limit N — > oo of the equations of Section 
4. To do this, we first define particle and hole densities as functions of the continuous 
variable x — 4*, 

p( x ) = ^ E - p h w = a7 £ *(* - $ ( 5 - 19 ) 

ie{i} m^{I} 

in such a way that p(x) + p h (x) — > 1 as iV — > oo. We can also write these in rapidity 
space by using the transformation rule for 5 functions, so that the Bethe equations 
become 

/oo 
d\'(f) 2 (\ - X')p(X') = 2nx(X) (5.20) 
-oo 

where we view x as an implicit function of A. Taking the derivative of this with respect 
to A and using ^jp- = p(A) + p h (\) yields 

/oo 
d\'a 2 (\- X')p(X') =p(A) +p h (X), AgR. (5.21) 
-oo 

For the particular case of the ground state, the occupation density pas(X) is non- 
vanishing in a symmetric interval [—B,B], with Pqs(X) vanishing. Outside of this 
interval, pes vanishes but not p^ s . A = ±B therefore represent the two Fermi points in 
the rapidity distribution of the ground state, which is obtained by solving 

PgsW + f B d\'a 2 (X - X')p G s(X') = ai (X), X G [-B, B\. (5.22) 

J-B 

The magnetic field dependence is encoded in the constraint 
B Ml 

JXp GS (X) = - = --a (5.23) 

where a is the field-dependent average magnetization per site along the z axis. These 
two equations determine B and pes-, and therefore also Pq S . We can write a formal 
solution as follows. Let us define the inverse operator L(X,X'), A, A' G [—B,B], inverse 
of the kernel in (I5.22p in the sense that 

-B 

dX'[5(X - A') + L(X, A')] [S(X' - A) + a 2 (A' - A)] = 5(X - A). (5.24) 

-B 

This operator is symmetric, L(A, A') = L(A',A), unique and analytic in its domain of 
definition [45]. In particular, the definition implies the identity 

a 2 (A- A) +L(A,A) + I dX'L(X, X')a 2 (X' — A) = 0, X,Xe[-B,B]. (5.25) 

J-B 
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In terms of this operator, we have the explicit solution of equation (15.221) for the ground 
state distribution, 

o (X) = l S -z dX ' [6{X ~ A ° + L(A ' A/)]ai(A ' } |A| " B ' (5 26) 

PGSy ' \ \\\>B. K ' ' 

Knowing p GS then yields p GS from (15.211) . namely 

P ^ (A) = { a,(A) - CdX'liX - \>) Pgs (\>) \x\ > B. ^ 

The ground state can also be obtained from the thermodynamic Bethe Ansatz formalism 
[45] in the following way. Given distributions p(A) and p h (X), the free energy / = 
(E — TS)/N is written to leading order in N as 

h r°° 

/=--+/ d\ [e p - T(p + p h )\n(p + p h )+Tp In p + Tp h In p h ] (5.28) 

in which we have suppressed the A functional arguments and defined the bare energy 

e (A) = /i-TrJsinCai(A). (5.29) 

Introducing the quasi-energy e(X) = T In , the condition of thermodynamic 
equilibrium SF = under the constraint of the Bethe equations (15.211) then gives after 
standard manipulations [13] (taking the limit T — > 0, so here and in what follows, e(X) 
is for the ground state configuration) 

e(A) + / dX'a 2 (X - A')e(A') = e (A), A e] - oo, oof. (5.30) 

J-B 

In particular, we have that 

e(±B)=0, e(A) < 0(> 0)forA e (£)[-B,B]. (5.31) 

Similarly to (I5.26p . we can also solve for e(A) = e~(X) +e + (X) with e ± (X) > (<)0 using 
the inverse integral kernel: 

I—B dX '[^ - A ') + |A| < B, 



e"(A) 



|A|>5, 



|A| < B, 

£ ( )_< > E (X)-f B B dX'a 2 (X-X')e-(X') \X\>B. (5 ' 33) 



The free energy simplifies to 



h 



B 



f = -- + J dA ai (A)e(A). (5.34) 
The magnetic equilibrium condition §f = then is 

By defining the dressed charge Z(X) as solution to 

Z{X)+ f dX'a 2 {X- A')Z(A') = 1, (5.36) 
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which we can solve as 

Z(X) = 1+ f d\'L(\, A'), (5.37) 

J-B 

we have the identity Z{X) = by making use of (15.291) and (15.301) . The Luttinger 
parameter K is given by the square of the dressed charge at the Fermi boundary (see 
e.g. [39]), 

K = Z 2 (—B). (5.38) 

The magnetic field dependence of the Fermi boundary B can be obtained by taking 
the h derivative of (15.301) : 

* dX'[5(X - X') + a 2 (X - A')]^ = J|. (5.39) 

Since e(—B) = 0, we have 

de{\) de{X) 
~~dX~ = ~ B ~ ~dB~\ x= ~ B 

and therefore 

dh_ _ e'(-B) 
dB ~ Z(-B)' 
The magnetization is 

df 1 f B ,..de(\) 1 rB 



(5.40) 
(5.41) 



a ' dh 2 



To get the susceptibility, we start from 

da_ f B ^ t ^dZ{X) 
dB ~ 



f dXai ^^T ~ a iW Z W ~ <~B)Z{-B). (5.43) 
The integral equation for the dressed charge (I5.36P gives 

^fP= - f B B dX>\5{X-X>) + L{\X')] 

x [Z(B)a 2 (X' -B) + Z(-B)a 2 {X' + B)} (5.44) 
which yields after simple manipulations and use of symmetry 

|| = -2p GS (-B)Z(-B). (5.45) 

The susceptibility is therefore given by 

da dBda _ PGS {-B)Z\-B) 

X dh dh dB e'(-B) ' 1 ; 
This expression will be related to the Fermi velocity after discussing elementary 
excitations (see equation (15.681) ). 

Finally, we will need the slope of the ground state rapidity distribution at the Fermi 
boundary, 8p g^ \-B- From the integral equation for pes, we can write 

= p GS (-B)[L(X, B) + L(A, -B)]. (5.47) 



Dynamical structure factor at small q for the XXZ spin- 1/2 chain 29 

This can be related to the derivative of the dressed charge by using the representation 

?§l = L{X ,B ) + H X,-B ) + J_y a -^l. (5,8) 

From the definition of L(X, A'), we can show that 
dL(X, X' 
dB 

and therefore 

dZ(X) 
dB " 

finally yielding 



= L(X, B)L(X', B) + L(X, -B)L(X', -B) (5.49) 
[L(X,B) + L(X,-B)]Z(-B), (5.50) 



dposW, _ p GS (-B) dZ(X) l 

dB " Z{-B) 8B l ~ B - (5 - 51) 

We will make use of these identities later, while relating the width of the two-particle 
continuum to field-dependent physical quantities. 

Let's now construct an excited state over the finite-field ground state by generating 
a single particle-hole pair. That is, we select a quantum number I p ^ {/ } associated 
to a particle and Ih € {I GS } associated to a hole, and write the excited state densities 
in x space as 

p(x) = pgs(x) + ±6(x - |) - ±6(x - §), 

P h (x) = PgsW ~ jj8(x - |) + ±5(x - |), (5.52) 

with once again p(x) + p h (x) — > 1 as N — > oo. We can again map to rapidity space, with 
X p < — B and \ Xh\ < B. Upon creating such a particle-hole pair, the induced distribution 
p(X) will be only very slightly shifted (order 1/N) as compared to the ground state one 
(for A ^ X p , Xh). We therefore define a backflow function K(X; X p , Xh) ~ 0(N°) as 

p(X) = Pgs(X) + 1 [K(X; X p , X h ) + S(X - X p ) - S(X - X h )} . (5.53) 

By subtracting the equations for the ground state from those of the excited state, the 
backflow function is shown to obey the constraint 

K{X; X p , X h )+ [ dX'a 2 (X-X')K(X'; X p , X h ) = -a 2 (A-A p )+a 2 (A-A ft )(5.54) 

J-B 

for A G [-B, B], with K — outside of this domain. We can again formally solve for K 
by applying the inverse integral operator 1 + L, 

K{X; Xp, Aft) = -a 2 (A - A p ) - / dX'L(X, X')a 2 (X' - X p ) - L(X, X h ). (5.55) 

J-B 

In terms of this kernel, the energy of the excited state is 

/oo pB 
dXe [p - Pgs] = e (Xp) - e (X h ) + / dXe (X)K(X; X p , X h ), (5.56) 
oo J-B 
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which can be rewritten after basic manipulations as (|A P | > B and \Xh\ < B) 

E-E GS = e(X p )-e(X h ). (5.57) 
Similarly, the momentum of the excited state is 



P-Pas= -MK) + M*h)~ f dXM^)K(^K^h)- (5.58) 

J-B 

Single particle-hole pairs as described above constitute a set of two-particle 
excitations labeled by the particle and hole rapidities X p and A^. This continuum is 
well-defined and spanned by the intervals X p G ] — oo, —B], Xh G [— B, B]. Assuming 
that the mapping from (X p , Xh) to (u, q) is one-to-one and onto and that the particle 
dispersion curvature is greater than the hole one (this monotonicity assumption will be 
discussed further in section 15731) . the highest energy stcite cit a given fixed momentum q 
will be given by the choice X p = X p (q), Xh = —B, where X p (q) is solution to 

q = -MK(q)) + M-b) + [ dtyi(A)#(A; x p ( q ) ; -b). (5.59) 



J-B 

Similarly, the lowest energy state will correspond to the choice X p = —B, Xh = Xh(q), 
where Xh(q) is solution to 

q = MMO)) ~ M-B) - [ dX(pi(X)K(X; -B; X h (q)). (5.60) 

J~B 

As discussed in Section 4, this continuum is well-defined (i.e. finite real solutions to 
both (15.591) and (15.601) can be found) as long as q < Mm(2k F , k^), with 2k F = 7r(l — 2a) 
and fcoo = 2a(n — (). This is illustrated in Figure (Q. The width of the two-particle 
continuum defined by these excitations will thus be given by the energy difference 
between these two limiting configurations, namely 

W(q) = e(X p (q)) + e(X h (q)) - 2s(-B) = e(X p (q)) + e(X h (q)) (5.61) 

where we have used s(±B) = 0. These functions are exact in the thermodynamic limit, 
in the sense that they allow at least in principle to obtain the exact function W(q) for 
the momentum region where these excitations are defined. These coupled equations 
unfortunately cannot be solved explicitly at nonzero magnetic field (where B is finite). 
We can however obtain analytical results in the small momentum limit, where these 
excitations always exist in a finite region at finite field. 

At small momentum, we can expand the width at fixed magnetic field as 

W = qW {1) + q 2 W (2) + 0(q 3 ) (5.62) 

with coefficients given explicitly by 

= ^(e(X p (q)) + e(X h (q)))\ q=0 , (5.63) 
= \^{< X M) + e(X h (q)))\ q=0 . (5.64) 
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Figure 6. Highest and lowest energy two-particle excited states at fixed momentum. 
The straight line represents the interval A e [-B, B] within which the ground-state 
rapidity pgs(A) is nonvanishing. \ p and A;,, respectively represent the positions of 
the particle and hole rapidities for the highest (top) and lowest (bottom) two-particle 
excited states at a fixed value of momentum. 



Let us treat the linear term first. Considering that (15.601) also defines a function q(\h) 
we can rewrite the hole contribution to the coefficient as 



de(X h ) I 



d 

Q- q e{K{q))\ q =o ^ 

9X h 



-13 



-13 



The denominator is obtained from (15.601) as 



dq 



2'Kai(\ h 



3 dXMX)- m - B '' Xh] 



(5.65) 



(5.66) 



-B 

where we have used (15.551) . the symmetry of L and partial integration. This contribution 
is by definition related to the field-dependent Fermi velocity, namely 







dq~ K """ iq 27T \^-b+ p G s(X) 
In particular, this allows us to relate the susceptibility to the Fermi velocity and the 
dressed charge using relation (I5.46p . 

Z 2 (-B) = ixvx- (5.68) 
For the particle contribution to the linear term, we find similarly that J^- = 
-2ixp h GS {X p ). Since lim A ^_ B - p GS {\) = ]im x ^- B + Pgs(X), we also have ^s(X p )\ q=0 = 

v since e is smooth around this point. Therefore, in the 



lim 



e-'(A) 



—v. 



(5.67) 



2n — A— B -fc gW 

momentum expansion (I5.62|) for the width, the linear term vanishes: 



= ^(e(X p (q)) + e(X h (q)))\ g=0 



0. 



(5.69) 



The width therefore depends at least quadratically on momentum. To compute the 
coefficient of the quadratic term, we first note that given a function X(q) and its inverse 
q(X), the chain rule allows us to write 

&_ 

dq 2 





dq 




-2 

( 


'd 2 e(X) 


d 2 q 


de(X) 


q=0 


dX 


-B. 




K OX 2 


-b dX 2 


-b dX 



. (5.70) 
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From (I5.59P and (I5.60p . we have that the particle and hole parts are related through 



dq(X p ) 



dq(\ h 



d\ p 
d 2 q{\) 



-B 



dX 2 p 



d 2 q(\ h 



-B 



9*1 



(5.71) 
(5.72) 



-B 



so using (I5.70p for X p (q) and \h(q), we obtain that the quadratic coefficient of the width 
can be simplified to 

2dq 2 

While this expression for the width is an end in itself, it is much more enlightening 
to relate it to more physical quantities by making use of the identities derived earlier. 



2 (s(X p (q)) + e(X h (q))) 





dq 




~ 2 d 2 e(X) 




q=0 


dX 


-B- 


dX 2 


-B 



(5.73) 



Starting from 
get 



d 2 e{\) 
d\ 2 



d 2 e(X) 



d e(X) 
dB 2 



OX 2 



-2ixp 



b and using (15.671) together with (I5.5ip and (I5.68p . we 

(5.74) 



as 



-B 



1 v dx 
2xdB 



Putting this in (I5.73P and making use of (15.411) . (15.671) and (I5.68P again, this finally gives 

3 dv 1 v dx 




v 

KX 



2dh 2 x dh 



(5.75) 



Since we have the identity K = Z 2 (—B) = nvx, this coincides with (I3.69p . It also 
reduces to the formula derived in [32] for A C 1 by linearizing the Bethe Ansatz 
equations. While our derivation was done for the anisotropic chain in the gapless regime, 
the same calculation can be performed for the isotropic antiferromagnet by simply using 
the appropriate scattering kernels in the Bethe equations. This result is however limited 
to chains with finite magnetization, in view of the fact that the region of validity of the 
excitations we have used to compute the width collapses to zero when the field vanishes. 



5. 3. Comparison with numerical form factors 

In order to compare the field theory results with the dynamical structure factor for finite 
chains, we first fix the parameters of the bosonic model introduced in section 13.11 We 
do that by calculating v(A,h) and K(A,h) = irv (A, /z)x(A, h) numerically using the 
Bethe Ansatz integral equations in the thermodynamic limit. 7/_ and r] + are obtained 
by linearizing the field dependence of v and K around some fixed ho and using (13.691) 
and (I3.70p . As examples, we consider three values of the anisotropy, A = 0.25, A = 0.75 
and the Heisenberg point A = 1, at a fixed magnetization per site a = —0.1. Tabled] 
lists the values of the important parameters (we set J = 1). Note that b is negative for 
a < (m > 0) because K decreases as we approach half-filling [2] . 

As mentioned in section 14.21 we are able to calculate the exact transition 

I I I 1 2 

probabilities F 2 (q,uj) = (0 a) for finite chains by means of the Algebraic Bethe 
Ansatz [T9J [201 121] ■ Figure [7] illustrates a typical result obtained for finite anisotropy 
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Table 1. Parameters for the low-energy effective model for A = 0.25, A = 0.75 and 
A = 1 and finite magnetic field ho (in all cases the magnetization per site is a = —0.1). 



A 


/'o 


V 


K 


a 


b 


n- 


V+ 


0.25 
0.75 
1 


-0.414 
-0.652 
-0.791 


1.087 
1.313 
1.399 


0.871 
0.699 
0.639 


0.306 
0.271 
0.256 


-0.050 
-0.145 
-0.188 


0.356 
0.409 
0.397 


-0.095 
-0.449 
-0.690 
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Figure 7. Numerical form factors squared (transition probabilities) for states with 
momentum q = 27r/25, for a chain with N = 200 sites, anisotropy A = 0.25, 
magnetization per site a = —0.1. The energies of the cigenstates are rescaled by 
the level spacing of the bosonic states predicted by field theory. The on-shell states 
are the ones at luN/2itv = qN/2-K = 8. 



and finite magnetic field. In contrast with the free fermion case, we observe two main 
differences when we turn on the fermion interaction A: First, the form factors for 
the two-particle (on-shell) states become uj- dependent; second, the form factors for 
multiparticle states are now finite and account for a finite spectral weight extending 
up to high energies. For the four-particle states (two particle- hole pairs), we expect 
(0|S'||a) ~ O(A), but this is not true near to fa vq where perturbation theory in 
the interaction diverges [28]. Figure [7] suggests that most of the exact form factors 
evolve smoothly from the XX point, except close to the lower and upper thresholds. If 
that is the case, the two-particle states still carry most of the spectral weight. In the 
thermodynamic limit, F 2 (q,u) has to be combined with the density of states factor 

D ^ = ^J2 6 ^- E * + E gs): (5-76) 
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Figure 8. Form factors squared for the two-particle states for two values of system 
size N (we set q = 27r/25, A = 0.25 and a = —0.1). The points seem to collapse 
on a single curve, showing very little size dependence. The minimum and maximum 
energies converge to the thresholds of the two-particle continuum when N — > oo. 



to define the lineshape of S zz (q, u>) (see ( 12. 3p ). 

We can count the states at each energy level of the finite system in the Bethe 
Ansatz the same way we count states for weakly interacting fermions. For example, 
in figure [7] we see n = qN/27r = 8 two-particle states with F 2 ~ 0(1). One can also 
verify that for n = 8 there are 14 states with two right-moving particle-hole pairs (of the 
form Cp l+ R c pi RCp 2+g2 R c p2 R \0)) and no states with three or more pairs. The 14 on-shell 
states with F 2 (q, uj) < 1CT 3 in figureHare all four-particle states. Furthermore, for small 
A the main contribution to the high-frequency tail (£ = uN/2nv > n) is due to states 
containing two particle-hole excitations created around the two different Fermi points 
[27] . If the momenta of the pairs at the right and left branches are q\ = 2irni/N > 
and q2 = 2irn 2 /N < 0, such that n\ = (£ + n)/2 and n 2 = —(£ — n)/2, then the 
number of such states is given by \n\ x n 2 \ = (£ 2 — n 2 )/A. This is in agreement with the 
counting of states in figure [7J We also find much smaller form factors for states with 
three particle-hole pairs (not shown in the figure). 

We now focus on the two-particle states inside the peak, with uo m vq. If we seek 
only these states with dominant form factors it is possible to reach much larger system 
sizes (we go up to 7000 sites). The number of two-particle states is always n = qN/2it. 
Figure [8] shows F 2 (q,u) for a fixed value of q = 2n/25 and two different system sizes. 

We extract 8u q from the numerical form factors as follows. We see from figure 
[8] that the separation between energy levels inside the peak is of order 8uj q /N and 
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Figure 9. Finite size scaling of the minimum energy for two-particle states with 
q = 2tt/25, A = 0.25 and a = -0.1. 

Table 2. Effective inverse mass, defined as the coefficient of the q 2 scaling of the width 
Su) g . The data are for a = —0.1 and anisotropy parameters A = 0.25, 0.75, 1. 



A 


l/m* FIT 


V- 


i (1+ 2A s infc F ) 

m V 7r r J 


0.25 


0.354 


0.356 


0.356 


0.75 


0.408 


0.409 


0.449 


1 


0.396 


0.397 


0.496 



decreases from ool{q) to uiu{q). As N increases, the maximum and minimum energies 
(N) converge to fixed values which we identify as the thresholds of the two- 
particle continuum. Figure [9] shows the finite size scaling of the minimum energy for 
A = 0.25, a = —0.1 and q = 2ir/25. The same A r_1 dependence is observed for the 
maximum energy. We use this scaling to determine the lower and upper thresholds 
u>L t u (q) in the thermodynamic limit for several values of q. We then calculate the 
width 5u q = u>u (?) — Wl (<?)■ As expected, we find that 8ui q = q 2 /m* for small q 
(figure [TO]) . Table [2] compares the coefficients \jm* FIT obtained by fitting the data with 
the predicted values of r\- taken from Table [U The perturbative result in (13.441) is 
also shown for comparison. The agreement supports our formula for the width in the 
strongly interacting (finite A) regime. Note that 7?_ is a nonmonotonic function of A. 

In figure[UJ we confirm that, despite the enhancement (suppression) near the lower 
(upper) threshold, F 2 (q,u) converges to the constant value F 2 (q,u) = K in the limit 
q — > 0, as expected from the box-like shape shown in figure [5] (see however the subtleties 
about the thermodynamic limit in section [872]) . This is in agreement with the fact that 




Figure 10. Width of the on-shell peak (based on the two-particle contribution) as a 
function of momentum q for a = —0.1 and two values of anisotropy: A = 0.25 (blue 
diamonds) and A = 0.75 (red circles). The lines are the best fit to the data. 



the exponents of the singularities at the edges are linear in q for h ^ [28] . Notice that 
F 2 (q,u) (and therefore S zz (q,oj)) is not a scaling function of (cu — vq) /8ui q . 

If the density of states D(q,u) for the two-particle states were constant, the 
extrapolation of figure [8] to the thermodynamic limit would be representative of the 
lineshape of S zz (q,u). This would be exactly the case if the exact energies E a — Eqs 
could be written as the sum of the energies of particles and holes with parabolic 
dispersion (as in a Galilean-invariant system, e.g. the Calogero-Sutherland model |29j). 
This is also the case considered in [28]. In our case D(q,u) does vary inside the peak 
because of the cubic terms in the dispersion of the particles in the Bethe Ansatz. For 
large enough N we can include the density of states factor ( 15.761) if we rescale F 2 (q,u) 
by the separation between energy levels inside the peak 

S»(q, u) = D(q, u)F\q, u) » ^ , (5-77) 

where Ej and Ej + i are the energies of the two-particles states, ordered in energy, with 
Ej — Eqs — The approximate density of states calculated this way is illustrated 
in figure [12l The resulting lineshape is shown in figure [131 This lineshape should 
be contrasted with the free fermion result in figure [H The exact boundaries of the 
two-particle continuum (dotted line in figure [13]) are actually shifted to lower energies 
relatively to the prediction ujj,L(q) = vq±r]^q 2 /2 (dashed lines) because of the cubic term 
in the exact dispersion, which was neglected in the field theory approach. Notice that 
there appears to be a peak at the exact lower threshold of the two-particle continuum. 
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(co-vqyScOq 

Figure 11. Frequency dependence of the form factors squared for N = 6000, A = 0.25, 
a = —0.1, and three values of momentum. The dashed line represents the field theory 
prediction F 2 (q,u>) = K ~ 0.871, as in figure [5] 




Figure 12. Density of states D(q,oj) for the two-particle states obtained using (|5.77l) . 
As in figure [Til we use N = 6000, A = 0.25 and a = —0.1. The main graph is for 
qN /2tt = 80. The inset shows the density of states for a smaller value of momentum, 
qN/2TT = 40. The solid lines are meant to illustrate the deviation of D(q,u>) from the 
linear dependence in uj. 
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Figure 13. Lineshape of S zz (q,u>) estimated from the two-particle states (Sf J (q,u>) 
in the notation of section l8?2|) . For this graph A = 0.25, a = —0.1, N = 6000 and 
qN/2-K = 80. The dashed line is the flat distribution of figure The dotted lines are 
the exact boundaries of the two-particle continuum in the thermodynamic limit . 



The result of Pustilnik et al. predicts that there is actually a power-law singularity at 
ul(q), which is related to the physics of the X-ray edge problem [2H]. We do not attempt 
to study the singularity in the form factors in this paper (see discussion in section . 
Interestingly, however, the density of states competes with the energy dependence of the 
form factors, leading to a minimum in S zz (q, u) above u>l{q) and a rounded peak below 
u>u(q). In the limit q <C cot hp we can linearize the density of states for the two-particle 
states 



D(q,u) 



2n/N 

E j+1 ~ E 'J 



1 



r]_q I 



1 + 



7g u — v q 



V- 



5u 



Q J 



(5.78) 



where 7 is a fitting parameter analogous to 7 in 02.21 j) for the free fermion model. The 
inset of figure fT2l shows the density of states for a smaller value of q = 2tt (40/6000). 
We have checked that D(q,u) becomes more linear and 7 converges to a finite value as 
q decreases. For A = 0.25 and a = —0.1 we estimate 7 ~ 1.11, which is larger than 
the value for free fermions 7 = sin(27r/5) ~ 0.951. Combining this density of states 
with the power-law singularity proposed in [28], the behavior near the lower threshold 
is described by the function 



S zz (q,co) 



K 

V-Q 



7g jqu-u L (q) 



2r/_ 



Suj n 



to 



Suj n 



(5.79) 



where \x q is the exponent of the X-ray edge singularity. The position of the minimum is 
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Figure 14. Lineshape for the Heisenberg chain at finite field (A = 1, a = —0.1, 
N = 6000 and qN/2-K = 80). Lines and symbols are represented as in figure fT3l 



then 



SUJ„ 



(5.80) 

for 7g/r/_ <C 1 and fi q <C 1. Since n q oc q for small q, the right-hand side of (15.801) 
becomes constant in the limit q — > 0. In this sense, the X-ray edge singularity and 
the energy dependence of the density of states are effects of the same order in q. We 
notice that the difference AS ZZ , defined between the maximum and the minimum of 
S zz (q,u), converges to a finite value as q — > (as it did for free fermions). The precise 
value depends on both the density of states D(q,u>) and the frequency dependence of 
F 2 (q,u) (which is approximately linear with a negative slope for \u — vq\ 5u q ). As 
a result, the slope of S zz (q, to) near the center of the peak diverges as 1/q 2 as q — > 0. 
This is a rather singular dependence of the lineshape on 7 and is potentially important 
for systems in which the dispersion is not exactly parabolic (e.g. due to band mixing in 
semiconductor quantum wires). 

Figure [14] shows the lineshape for the isotropic point A = 1 and the same values of 
a and q used in figure [131 In comparison with the weak coupling value A = 0.25, there 
is an enhancement of the singularities near the lower and upper thresholds. The shift of 
the peak to lower energies (another "g 3 effect" ) is also more pronounced, but the width 
is very well described by the field theory formula (prefactor given in table [2]). 

Finally, let us comment on the validity of the q 2 scaling for the width as a function 
of magnetic field. Figure [15] shows the dependence of the coupling constants of the 
irrelevant operators on the magnetization a for A = 0.25 and A = 1. From the field 
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Figure 15. Parameters ij± for the low energy effective Hamiltonian as a function of 
the magnetization a < for two values of anisotropy: (a) A = 0.25; (b) A = 1. For 
a > 0, we have r?±(— c) = — T)±(cr). 



theory standpoint, we expect that the q 2 scaling is valid as long as r]_q 2 <C vq (the peak 
is narrow) and t]±q 2 ^> ^q 3 (the cubic terms yield the leading correction to the free 
boson result and we can drop operators with dimension four and higher in the effective 
Hamiltonian). For A = 0.25, we see that rj± follow the behavior predicted by the weak 
coupling expressions (13.441) and (13.451) . vanishing at o = 0. In this case, r)± are always 
of 0(1). The restrictions are similar to the ones for the approximation (I2.24p for the 
dynamical structure factor of the XX model, namely q <C &p and q <C cot hp (which 
becomes q <C ira for small a). On the other hand, for A = 1 we find that \rj±\ — > oo 
as o~ — > 0. This is a direct consequence of formulas (13.691) and ( 13.701) in the strong 
coupling regime. It is known that the magnetic susceptibility at small fields is given 
by x(h) ~ const + C\h 2 + C^/i 8 ^ -4 , where are constants [43l I4H] . The exponent 
8K — 4 is a manifestation of the Umklapp scattering term at zero magnetic field. As a 
result, dx/dh diverges as h — > for K < 5/8 or A > cos(7r/5) ~ 0.81. In other words, 
the Luttinger parameter has an infinite slope at h — (see [37] for the isotropic case). 
Since dK/dh and dv/dh have opposite signs, t]_ goes through zero for a finite value of 
a. Therefore, we predict that Su q is a nonmonotonic function of a for A > cos(7r/5) and 
| a | 1. The sign change in 7/_ is reflected in the Bethe Ansatz data as the inversion of 
the ordering of the energies of the two-particle states as a function of hole momentum. 
For A = 1, the "inversion point" where 7/_ = occurs at \<Ji nv \ ~ 0.030. At this point 
we observe that Su q oc q 3 for q na. However, r) + ~ 0(1), so the lineshape defined 
by the two-particle states must be different from the one at zero field, where there is 
also a q 3 scaling (see section [7J. We have also confirmed that the q 2 scaling holds in 
the region where < and q <^ n\a\ (figure [T6|) . In this regime we find that the sign 
change of 7/_ is accompanied by the inversion of the lineshape of S zz (q,u>): The form 




Figure 16. Width 5ui q for A = 1 and a = —0.01 (in the region where rj- < 0). The 
blue diamonds represent the width defined as the difference between the maximum 
and minimum energies of the two-particle states calculated in the Bethc Ansatz, 
extrapolated to the thermodynamic limit. For q -C 7r|tr| , we recover the behavior 
5ujq = \r]-\q 2 , with \r]-\ w 0.246 (dashed line). 



factors appear to vanish at the lower threshold and are peaked near the upper threshold 
(with a possible divergence at ujj) (figure [TTj) . There is no maximum or minimum near 
the edges in this case. In order to understand this result, we recall that a converging 
X-ray edge is possible in strongly interacting systems. An important point is that the 
exponents of the X-ray edge singularities calculated in [28], which predict a diverging 
X-ray edge, are valid only to first order in the interactions. Second order corrections, 
which usually have the opposite sign because of the orthogonality catastrophe, tend to 
kill the singularity at the lower edge [48]. For even smaller values of a, the divergence 
of rj- seems to be consistent with the Bethe Ansatz results. For \a\ < <7j nw , the width 
increases as \a\ decreases at least down to a = —0.001, the lowest magnetization we 
were able to analyze. In the limit a — > and \r]±\ ^> 1, we expect for the isotropic point 
(using the results of [4^||47]) 

— -> V- - 7= J , (5-81) 

3 8v2crln \(Jo/a\ 

where do = ^/32/7re. According to the conditions ^q 2 <C r]_q ^ v, the field theory result 
which predicts the q 2 scaling for a small fixed q breaks down both near the inversion 
point <Ji nv and for a — > 0. In the limit a — ► 0, as mentioned in section 4.2, the set of 
allowable quantum numbers for the single particle-hole excitations becomes empty, as 
the Too quantum number tends to N/2 (meaning that the particle part becomes trapped 




Figure 17. Lineshape of S zz (q 1 Lo) for magnetization below the inversion point, i.e. 
\a\ < <r mv (A = 1, a" = -0.01, N = 7000 and qN/2n = 14). This value of q is in the 
domain where 5ui q ~ g 2 (see figure IT5|) . but the lineshape is inverted. In contrast with 
figurc[8j the form factors (shown in the inset) are peaked at the upper threshold of the 
two-particle continuum. 



at the Fermi surface), and this family of excitations disappears. The vanishing field two- 
particle continuum at nonvanishing momentum is then obtained from considering the 
next simplest excitations, which are states having two holes (spinons) within the ground 
state configuration together with a single negative parity one-string (or, for the XXX 
chain, an infinite rapidity). At finite but small field, the contributions from these states 
dominates S zz (q,u) for q ^> na and allows to smoothly recover the zero field behavior. 
A full discussion of all the possible lineshapes as a function of A and a together with 
the characterization of the dominant families of excitations is accessible from the results 
of [21], but is beyond the scope of the present paper. 

6. High-frequency tail 

We now turn to the calculation of S zz (q, lj) in the frequency range 7 9 <C oj — vq <C J, 
where finite order perturbation theory is expected to be valid. This off-shell spectral 
weight is possible because the rj + interaction allows for two-boson intermediate states 
with total momentum q = qi + q2 but energy oj = v \q\\ + v \q2\ > v\q\ if sign(gi) = 
— sign(g2)- In other words, the incoming boson can decay into one right-moving and 
one left-moving boson, which together can carry small momentum but high energy 
to 3> v \ q\. In the limit A < 1, this is equivalent to a state with two particle-hole pairs 
created around the two different Fermi points [27]. In this sense, our t] + is analogous 
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Figure 18. Diagrams at 0{rn_) for the calculation of the tail. 

to the U q interaction in [25] • We should stress that, although the tail carries a small 
fraction of the spectral weight of S zz (q,uj), it is important for response functions that 
depend on the overlap of two spectral functions, e.g. the drag resistivity in the fermionic 
version of the problem [27|. In our formalism the calculation of the tail provides a direct 
quantitative check of the accuracy of the low energy effective model against the form 
factors calculated by Bethe Ansatz. 



6.1. Field theory prediction 



The lowest-order correction to x (9> ^n) due to the r] + interaction is 



Sx {q,iuJ n ) 



dxe~ 



iqx 



dre^S X (x,r), 



(6.1) 



where 5x {x, r) is the correlation function in real space given by 

2 



S X {x,t) = — ~ 

71 Z 



6 



d X\ \ d %2 



x {d x <j> (x) [{d x ip L (l)) 2 d x <p R (1) -(R^L) 
x [(d x <p R (2)) 2 d x <p L (2) L)] d x <f> (0) .) 

This corresponds to the diagrams in figure [18j 5x can be factored in the form 



Sx (q,iu) 



K 
27 



D%\q,iu)+D™(q,iu) U RL (q, iu) 



(6.2) 
(6.3) 



where H R l (q, iuj) is the bubble with right- and left-moving bosons 
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After integrating over frequency, we get 

k(q + k) 



(6.4) 
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(6.5) 
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where A ~ kp is a momentum cutoff. Note that the real part of IT^ is ultraviolet- 
divergent, but the imaginary part is not. The integration over the internal momentum 
yields 



, ril j A 2 uj 2 - v 2 q 2 Uygf-^ + i^f 



S8»{ q , U ) = ^- v -^^L. (6.7) 



9 [2t; 8t> 3 ( ^ 6 ' 6 ' ) 

Finally, using Eqs. (I6.3P and (12. 4p . we find that the high-frequency tail of S zz (q, uj) 
is given by 

Kr] 2 + q A 9{uj - vq) 
18u uj 2 — v 2 q 2 

This is the same uj~ 2 dependence obtained for weakly interacting fermions with parabolic 
dispersion [27] . Since the small parameter is rj + ~ A/m, we approach the perturbative 
regime either by A — > or m — > 00 (more precisely, q/mv — ► 0). In this limit, our result 
(16. 7p agrees with equation (19) of [27] if we use (I3.45P and U (q) = (A/2) cosg. 

Since our model predicts that 5S ZZ (q, uj ^> vq) ~ O^r] 2 ^), one interesting 
consequence is that there will be no tail in S (q, uj) for models where the Luttinger 
parameter K is independent of particle density, since then 77+ = according to ( I3.70p . 
This is the case for the Calogero-Sutherland model, where K is a function of the 
amplitude of the long-range interaction only [50J. 

The divergence of the high-frequency tail of 5S ZZ (q, uj) as uj — > vq confirms that 
the on-shell region is not accessible by our standard perturbation theory in the band 
curvature terms. The matching of the tail to the on-shell peak at uJuil) is a complicated 
problem that has only been addressed in the regime A 1 (see [28]). The (uj — vq)" 1 
divergence in (16.71) comes from the frequency dependence of the external legs in the 
diagrams of figure [181 It is easy to see that if the bosonic propagators are replaced 
by the "dressed" propagator (all orders in t]_) given by (]5.12|) . the singularity at the 
upper threshold uJu(q) becomes only logarithmic. This supports the picture that the 
r] + interaction only modifies the shape of the on-shell peak very close to the edges. We 
expect that 77+ will contribute to the exponent of the singularity at the edges, since 
the exponent \x q derived in [28] picks up corrections of second order in the interaction 
between right and left movers, i.e. 0(r}+). As discussed in section [57T1 we believe that 
77+ does not affect the width to 0(q 2 ). Evidence for that is that the perturbation theory 
in i] + (second order given by (16.61) ) does not generate terms with the same q and uj 
dependence as in ( 15.91) and ( 15.101) . If the frequency dependence is regularized in the 
peak region by summing the perturbation theory in the diagrams involving rj + are 
always suppressed by higher powers of q because of simple kinematics. Inside the peak 
the energy of the left moving boson that is put on shell when taking the imaginary 
part of x(q,u) (as in the "unitarity condition" method used in [35]) has to be of order 
Sujg = n^q 2 or smaller, which constrains the phase space for the internal momenta. 
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Figure 19. Tail of S zz (q,u>) for q = 27r/50 and Suj q < w - vq < J. The red dots 
represent the sum of the numerical F 2 (q, u>) identified with each energy level predicted 
by field theory (c.f. figure [7|). The solid line is the field theory result (|6.10p . The chain 
length is N = 600. (a) a = -0.1, A = 0.25; (b) a = -0.1, A = 0.75. 



6.2. Comparison with numerical form factors 

For a finite system with size N, the result for 5S ZZ (q, u) must be expressed in terms 
of the transition probabilities F 2 (q,u). If the intermediate bosons carry momenta 
qi t 2 = 2imi t 2/N, such that qi + q 2 = q = 2irn/N, the energy levels are given by the sum 
of their individual energies u = v \qi \ + v \q 2 \, i.e., 
2irv£ 

u t =— r , £ = n + 2,n + 4,... . (6.8) 
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Thus field theory predicts a uniform level spacing Attv/N above the mass shell. It is 
easy to verify (by simply replacing the integrals by sums in momentum space) that 
SS ZZ (q, u) for the finite system can be written as 



where F 2 (q, uo) = | ^0 \S* | Q;)| 2 , with \a) a two-boson intermediate state, is the transition 
probability for the state with energy ue and is given by 

F 2 (q, u) = 2v 5S" (q, u>) = (6.10) 

We compare our field theory prediction with the form factors calculated numerically 
for a chain with N = 600 sites. We take q = 27r/50 (n = 12) and the previous values 
a = —0.1 and A = 0.25 or A = 0.75 (for which the parameters are shown in tabled]). As 
we saw in figure [7J the energies of the eigenstates calculated by BA are actually scattered 
around the values of uji predicted in (16.8)1 . The broadening becomes comparable with 
the level spacing 4irv/N when I « 30 (ou ~ 0.4J). Again the number of states agrees 
with a picture of multiple particle-hole excitations based on perturbation theory in the 
interaction. These features are not predicted by the bosonization approach. In order 
to make the comparison with (16.101) . we group the form factors that can be identified 
with a given energy level and plot the total F 2 (q, u>) as a function of the integers 
I = LdN/2irv. We emphasize that for very large £ we expect deviations from the lowest- 
order field theory result due to the effect of more irrelevant operators we have neglected. 
The results are shown in figure [19j 

7. The zero field case 

So far we have focused on the dynamical structure factor at finite magnetic field, which 
is somewhat analogous to interacting fermions with parabolic dispersion. One may then 
ask whether the field theory calculations can be applied to the case h = (kp = 7r/2). 
Let us first review what is known for the free fermion point A = 0. In this case S zz (q, u) 
is still given by ( 12.191) . but the thresholds of the two-particle (two-spinon in the Bethe 
Ansatz solution) continuum are given by 

(q) = J sing, (7.1) 
uu{q) =2Jsin|. (7.2) 

As a result, S zz (q, lj) develops a square root divergence at the upper threshold uju(q). 
The width now scales like q 3 for small q 

<H « If. (7-3) 

A crossover from q 2 to g 3 is observed as we decrease the magnetic field (or, equivalently, 
increase q hp) so as to violate (12.151) or (14. 9p . The result (j7.3p is also obtained by 
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keeping the leading correction to the linear dispersion around kp 

7/c 3 



R,L 



± vpk — 



6 



(7.4) 



where 7 = Vp — J. Bosonizing the band curvature term according to (13.37j) . we find 

24 



sn bc 

which can be rewritten as 

sn bc 



(dlfa) 2 :+{R^L), 



(7.5) 



-^:(^ R ) 2 ::(^ R ) 2 :+(i?^L), 



(7.6) 



as follows from the operator product expansion of (17.61) . 

In the interacting case we also have to keep track of the irrelevant interaction terms, 
including the Umklapp interaction in (13 .4p . The general form for the leading irrelevant 
operators for zero field is [JH E2] 



5H 



12 



I 2 : : (d x ip R ) 2 



+ 



■ {d x (p R ) 2 : : (d x <f L y 



(d x ip L ) 



: : (d x tp L ) ■] 



2ir 



cos 



(4VvrK0) + . 



(7.7) 



where the dots stand for higher dimensional local counterterms. The coupling constants 
to first order in A can be obtained from the bosonization of the band curvature term 
and the irrelevant interaction terms. We find 
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The exact coupling constants for finite A can be taken from 

v r (life) r3 (life) 

2K 



c 

Ai 



AttK r ( 

V 

tan 

2tt 



-M r 3 ( 

7VK 



y 



2K 



AvT{2K) 
T(l - 2K) 



£(j+^h) 

2 v^r (1 + ^L.) 



4X-2 



(7.8) 

(7.9) 
(7.10) 

(7.11) 



where v and K are given by (13.161) and (I3.17p . respectively. 

One important point is that the other possible type of dimension-four operator 
{dx'PR) 3 d x ipi + R <-» L is absent from the effective Hamiltonian for the XXZ model. 
We see this directly when calculating the coupling constants to first order in A, but 
we can also show that it remains true for finite A by imposing the constraint that the 
XXZ model is integrable [32]. Integrability implies the existence of nontrivial conserved 
quantities, the simplest one of which is the energy current operator J E = . j? given 
by J5HE2] 



j2 \ ^ I" qU C 2 Q x Q x Q z CV 1 A I Q x QV Q z Q z CU Q a 

j 

+A(S*_ 1 S*S« +1 - SV^Sj+J] . (7.12) 
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The latter is defined by the continuity equation of the energy density at zero field 

jf +l -jf = -d t H 3 = i[H^Hl (7.13) 

where H = . Tij is the Hamiltonian (12.11) with h = 0. One can then verify that J E is 
conserved in the sense that [J E , H] = 0. 

Let us now look at the corresponding quantity in the low energy effective model. 
In the general case, we consider the Hamiltonian density 7i = Hll + STi. + S7i 3 , where 
SH is given by (17. 71) and we also add the interaction 

5H 3 = vrC 3 [(d x (p R f d x (p L + (d x ip£f d x (p R ] . (7.14) 

We obtain the energy current operator from the continuity equation in the continuum 
limit 

d x f(x) = -d t H(x) =i J dy [H{x),H{y)\ . (7.15) 

The energy current operator for the Luttinger model (with (± i3 = Ai = 0) takes the 
form 

Jo = J dxj^{x) = V — J ' dx [(d x ip R f - (d x ip L ) 2 ] 

= -v 2 I dxd x <pd x 0. (7.16) 



This coincides with the spatial translation operator of the Gaussian model [53]. A 
nontrivial consequence of the conservation law arises when we consider the corrections 
to J E due to the irrelevant operators. We keep corrections up to operators of dimension 
four. Using (I715jl . we find J E = J E + 5J E with [32] 

SJ E = nv j dx jy [{d x ip R f - (d x ip L ) A ] 

+2C 3 [(d x ip R f d xVL - (d x ip L f d x ip R ] } . (7.17) 

Note that there are no first-order corrections to J E associated with the ( + interaction 
or the Umklapp scattering Ai. (The case of the Umklapp perturbation was discussed in 
[54]). The conservation of J E up to dimension-four operators implies 

[J E , H] = [J E , H LL ] + [J E , 5H] + [J E , 6H 3 ] + [5 J E , H LL ] = 0. (7.18) 

Since J E is conserved in the Luttinger model, we have [J E , Hll] = 0. In fact, J E 
commutes with any local operator of the form j dxO(x) under periodic boundary 
conditions [53]. As a result, [J E ,5H] = [ J E , SH^] = as well. We are left 
with the condition that the commutator [5 J E , Hll] vanishes. This is automatically 
satisfied by the contribution from the term because it does not mix R and L and 
[{d x f R ) 4 — {d x ifL) 4 , Hll] is a total derivative. We then have 

[5 J E , H LL ] = kv 2 ( 3 J dx J dx' x 

[(d x (p R ) 3 d x (p L - (d x ip L ) 3 d x <f R , {d x i(p R f + (d x xp L ) 2 ] 

= Amv 2 C, 3 j dx [{d x tp R f d 2 x ip L + {d x ^ L f d 2 x ip R ] . (7.19) 



Dynamical structure factor at small q for the XXZ spin- 1/2 chain 



49 



Therefore, [J E , H] = <£> Cs = 0. 

This argument also applies to the finite field case. The model is still integrable for 
h =/ 0. Although the relevant quantity for thermal transport is now a linear combination 
of the energy current and the spin current operator (which is not conserved for the XXZ 
model), the energy current operator given by (17.121) commutes with the Hamiltonian 
(12. ip for all values of h [52j[55j. The corresponding conserved quantity in the low energy 
theory is the current operator J E obtained from the effective Hamiltonian at zero field, 
which has no dependence on the coupling constants r)±. Clearly, [jQ,5H(h =/ 0)] = 
for SH(h 7^ 0) given by (13.431) . so integrability poses no constraints on the coupling 
constants r)±. 

We have checked that £3 7^ for a nonintegrable model obtained by adding to the 
XXZ model the following next-nearest neighbour interaction 

6H nnn = JA'J2SjSj+2> ( 7 - 2 °) 

3 

which is mapped by bosonization onto 



SH nnn = J A' / dx 



■!N 2+ y(^ 4 



(7.21) 



The first term in ( 17.2ip is quadratic in the bosons and modifies the velocity and the 
Luttinger parameter of the Luttinger model. The second term is the irrelevant operator. 
To first order in A and A', we find that it gives rise to a (3 term in the Hamiltonian, 
which is given by 

29 

5H nnn ~ — -J A' [(d x <p R f d x cp L + (d x ip L f d x ip R ] . (7.22) 

This shows that, unlike the XXZ model, a low energy effective model describing a 
nonintegrable model must in general contain the £3 interaction. 

This result establishes a connection between integrability and the field theory 
approach, by means of a restriction on the coupling constant of a band curvature type 
operator in the low energy effective model. More generally, if we keep more irrelevant 
operators in the effective Hamiltonian, integrability should manifest itself as a fine tuning 
of the coupling constants and the absence of certain perturbations. This connection may 
be important for understanding the role of integrability in the transport properties of 
one-dimensional systems [5Tj . 

With (3 = 0, only ( + and Ai mix right and left movers. We can apply second order 
perturbation theory in these interactions to calculate two contributions to the high- 
frequency tail in the frequency range Su q <w-i;g€ J |32j. For a finite chain with N 
sites and fixed momentum q = 2nn/N, the ( + operator gives rise to intermediate states 
with discrete energies uj^ = 2ttv£/N, £ — n + 2, n + 4, . . .. In the thermodynamic limit, 
the contribution to the tail is 

K{C + /vf 2 f^-yV 
192v q V V 2 

The states generated by the Umklapp operator have energies ug = 2ttv(£ + 4K)/N, 
£ = n, n + 2, ... . For 4itv/N ^> 5u q it is easy to separate this contribution from 
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the ( + one because of the shift in the energy levels by the noninteger factor AK. The 
corresponding contribution to the tail is 

") = Y^) {2v)3 ~ 8Kq2 ^ 2 - v2q2 ^ K ~ 3 9{uj - vq) - (7 - 24) 

The derivation of equations (17.231) and (17.241) . as well as the result for the finite system, 
is presented in the appendix. 

The result in ( I7.23|) and (17.241) shows that the Umklapp operator (dimension 
4K) yields the dominant contribution to the tail near u ~ vq. For < A < 1/2 
(3/4 < K < 1), the next-leading contribution is given by (dimension 4). For 
1/2<A< 1 (1/2 < K < 3/4) it is important to include in the effective Hamiltonian 
the operator 

5Hx 2 = \ 2 d x 6cos(4V^K(f)) J (7.25) 

which is a descendant of the Umklapp operator and is allowed by all symmetries. This 
operator has dimension 4K + 1 and is less irrelevant than ( + for K < 3/4. Another 
reason to include A 2 is that for K < 3/4 the exact amplitude ( + in (I7.10p diverges at the 
points K = 1/2 + l/(4n), n > 1. This divergence has been discussed in the context of 
corrections to the bulk and boundary susceptibility of the open XXZ chain [46J . There 
it was found that the susceptibility as a function of magnetic field or temperature has 
a correction of first order in ( + . However, the corrections at any given order of h or 
T are always finite because the divergences of ( + at the points K = 1/2 + l/(4n) are 
cancelled by the contribution from the Umklapp operator and the cancellation gives 
rise to logarithm corrections. In our case the tail (I7.23P is of order so it must be 
cancelled by a more irrelevant operator with a independent dimension. Note also that 
the 0(\1) term in (I7.24p does not diverge. Therefore, in order to recover a finite high- 
frequency tail for all values of < A < 1 it is necessary that the amplitude of the terms 
generated by A2 also diverge (have poles) at the above values of K. By simple power 
counting, we expect that the divergence of the £+ term at K = 1/2 + l/(4n) will be 
cancelled by the term of second order in A2 and 2(n — l)-th order in Ai, which scales 
like 5S ZZ ~ \\ {n ~ l) \ 2 2 q 2 {uj 2 - uV) n(4K_2) - For n(AK - 2) = 1, this term has the same 
q and u dependence as in (17. 23ft and the cancellation is thus possible. In principle it is 
possible to determine the amplitude A2 from the Bethe Ansatz, but that would require 
solving the Wiener- Hopf equations to higher orders than was done in [4~6] . 

Computing the broadening 5u q for h = from bosonization is much more 
challenging. Since is the only vertex present at the free fermion point, the naive 
expectation is that we could derive the renormalization of the width at zero field by 
summing all orders of as we did for rj_ in section 15.11 Although we now have to 
deal with a four-legged vertex, which introduces three-boson intermediate states, the 
calculation of the lowest order diagrams is not much harder than the finite field case. 
However, the fundamental difference is that for h = the broadening has to be produced 
by dimension-four operators and is therefore of the same order of q as the changes in the 
lineshape (i.e., the density of states factor and the singularities near the thresholds). 
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Figure 20. Numerical form factors squared for A = 0.25 and zero field. The chain 
length is N — 600 and the momentum is set to q = 27r/50. 



That implies that the lineshape of S zz (q, uj) for A ^ cannot be approximated by the 
free fermion result in the limit q — > 0. Therefore it is not clear what the expansion of 
bosonic diagrams should sum up to. 

Figure [201 shows the form factors for A = 0.25 and h = calculated numerically by 
Bethe Ansatz for a chain of N = 600 sites. In agreement with the field theory prediction, 
the states in the high-frequency tail cluster around the energy levels ujN/2ttv = i 
(corresponding to the ( + contribution) and ujN/2iiv = £ + 4K rs £ + 3.45 (the Ai 
contribution). The comparison between the Bethe Ansatz data and the field theory 
results for the tail at zero field was done in [32], confirming the validity of (17.231) and 
(I7.24p for b~oj q < w - «g < J. An important consequence of integrability for the 
lineshape at zero field is that for A < 1/2 (K > 3/4) the tail decreases as uj — * vq. A 
finite ^3 interaction would produce a contribution to the tail that diverges as uj — > vq, 
similarly to the finite field case (see appendix). There is no such contribution in the 
Bethe Ansatz data for small A. 

To study the peak region we can focus on the two-spinon states only (with uj ps vq 
and form factors of 0(1)) and reach lengths up to N = 4000. We see that F 2 (q,u) 
is dominated by the two-spinon contribution, except very close to the upper threshold, 
where that contribution vanishes (inset of figure l2TT) . Unlike the finite field case, the 
rescaled F 2 (q,u>) does not become flat in the limit q — > 0. The density of states for the 
two-spinon states is known exactly 

1 
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qN/2% = 40 




[co-co(q)]/8co 

Figure 21. Two-spinon contribution to S zz (q,oj) at zero field for N — 4000, 
A = 0.25 and three values of q. The number of two-spinon states is given by 
qN/An. Inset: transition probabilities F 2 (q,uj) used to calculate S zz (q,u>). We denote 
5u q = uJu(q) -Uh{q) and Q(q) = [uu(q) +u; L (q)]/2. 

where ouu(o) — 2vsin(g/2). The two-spinon contribution to S zz (q,u) for zero field and 
A = 0.25 obtained by multiplying F 2 (q,u) by the above density of states is shown in 
figure EU 

We know from the exact solution for the two-spinon dynamical structure factor at 
the Heisenberg point A = 1 that there is a square-root divergence (with a logarithmic 
correction) at ul{q) an d that the same contribution vanishes at uijj(q) [56] . Such 
behavior is completely opposite to what happens at the free fermion point (see [6]). 
Numerical results suggest that the exponents change smoothly from A = to A = 1, 
with spectral weight being transferred from the upper threshold to the lower threshold 
as A increases [2HIMI25] . 

In addition, the renormalization of 5u q defined as the width of the two-particle 
continuum in the Bethe Ansatz solution is known exactly. The thresholds of the two- 
spinon continuum for < A < 1 are a simple generalization of ( 17. II) and ( 17.2ft . with J 
replaced by the renormalized velocity v given by (I3.16P [39l [23] . As a result, the width 
for finite A is 




(7.27) 



Since v ~ J(l + 2A/n) for A C 1, the above expression is different from the 
renormalization of Therefore, the renormalization of 5u q by interactions is not 
given by 
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A proper treatment of the dimension-four operators which allows to predict the 
lineshape of S zz (q, u) at zero field remains an open question. 

8. Sum rules and finite size effects 

8.1. Sum rules 

We checked the accuracy of the Bethe Ansatz data by calculating the following sum 
rules 

Ho) = J^S^u) = ^{S*Sl q ) (8.1) 

and 

L (q) = r — uS zz (q, u) = -2^- sin 2 1 . (8.2) 
W J 2tt W ' ; N 2 V 1 

These sum rules can be expressed in terms of sums over the form factors calculated by 
Bethe Ansatz for finite chains as 

w?) = ^EK°l^l a >l 2 ' ( 8 - 3 ) 

a 

LBA(q) = ^ ^ ~ Egs ^ K° 1^1 a )\ 2 ■ ( 8 - 4 ) 

a 

The identity in (18.21) is a consequence of the /-sum rule. The first moment sum rule 
L(q) can then be calculated exactly by using the BA result for (H xy )/N = 2(SjSj +1 ) = 
eo — deo/dA, where eo is the ground state energy per site [57] . Since there are no exact 
results for I(q), we first compare Iba with the lowest-order field theory result in (15.161) 

I(q) « I FT {q) = (8.5) 

This should be a reasonably good approximation for small q. 

The other possibility is to calculate the static correlation function by DMRG. In 
Table [3] we show DMRG results for periodic XXZ chains. The results used the standard 
DMRG finite system method [581 ES, EQ], but with extra noise terms added to the 
density matrix to speed convergence in the number of sweeps for the more difficult 
periodic boundaries case [61]. We see that for these measurements the finite size effects 
are very small for N = 100, and that the truncation error depends significantly on N. 
We can obtain results for I(q) to an accuracy of 10 -6 or 10~ 7 by using the m = 2400 
results for N = 100. Finite size corrections for larger N appear to be roughly the same 
size. 

A comparison between the sum rules obtained for the BA data and the values 
expected from the equations above is shown in tables H] and [5j In all the cases shown 
here the Bethe Ansatz agrees with the DMRG and the exact results to better than 0.1%. 
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Table 3. DMRG results for I(q), for A = 0.25 and zero field (a = 0). The truncation 
error is e and m is the number of states kept per block. Between 10 (m = 1200) and 
14 (m = 2400) sweeps were performed. 



N 


/(27r/50) 


/(27r/25) 


e 




ID 


50 


0.017237038 


0.034518387 


7.5 x 10" 


13 


1200 


100 


0.017237380 


0.034518764 


9.4 x 10" 


11 


1200 


100 


0.017237138 


0.034518522 


2.1 x 10" 


-11 


1600 


100 


0.017237109 


0.034518505 


1.5 x 10" 


-12 


2400 


200 


0.017237364 


0.034518657 


1.2 x 10" 


-10 


2400 



Table 4. Sum rules for A = 0.25 and zero field (a = 0). First sum rule: Results for 
the BA data, Iba, with a chain length of N = 400 compared with the field theory 
approximation (|8.5p . Ift, and results from DMRG, Idmrg, f° r a chain with N = 100 
sites. Second sum rule: Results for the same BA data, Lb a, compared to the exact 



result ((57 


^ : -^cxact ■ 










Ift 


Idmrg Lba Lexact 


y 50 

a - — 

y 25 


0.017236 

0.034513 


0.017228 
0.034457 


0.017237 0.002498 0.002498 
0.034518 0.009952 0.009953 


Table 5. 


Same sum 


rules as i 


in table HI but for finite field (a = - 




Iba 


Ift 


Idmrg Lba Lexact 


y 50 

a= 32 
y 25 


0.017418 
0.034880 


0.0174 
0.0348 


0.017419 0.002376 0.002378 
0.034883 0.009468 0.009474 



8.2. Size dependence of the n-particle contributions to S(q,u>) within Bethe ansatz 

For any finite N, S(q,u) is a sum of 5-functions. However, for N 3> 1, we expect to be 
able to approximate it by a continuous function of uj. One way this can be done is by 
"binning"; i.e. we can define: 

%,a;)^^^'|<0|^|a>| 2 , (8.6) 

a 

where the sum is restricted to energies such that 

uj- Alo/2 < E a <uj + Alo/2. (8.7) 

The bin size is chosen so that 2%/N Auj/v 2tt. For a fixed N, S is defined in ( 18.61) 
for discrete values of q. However, an extra binning could be defined for the wavectors 
as well, so that q could be held fixed as we increase N. Alternatively, we can use (18.61 ) 
and increase by an integer factor to reach the limit N — ► oo in such a way that the 
ratio qN/2n is always integer. In the large A^ limit S(q,u>) becomes a smooth nonzero 
function normalized so that: 

1 v-^ f du „, . 1 . 
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We can define an approximation to this function, for large but finite N, S^(q, u), using a 
fixed (small) bin size. Furthermore, we can decompose this function into the 2n-particle 
contributions: [t| 



It is convenient to discuss the finite size behavior in terms of these binned functions. 
While Sj^(q,u) must have a finite limit as N — > oo, it is possible that each single 
contribution «S|p(g, uj) vanishes. The question whether or not this is the case has 
been addressed for the isotropic antiferromagnet in zero magnetic field. Here it is 
known that two-spinon [TBI [15] as well as the four-spinon [18] contributions are finite 
in the thermodynamic limit. For the A = case, on the other hand, the two-spinon 
contribution to the transverse structure factor S + ~(q, u>) vanishes in the thermodynamic 
limit |56j . 

For the XXZ chain in a finite magnetic field the form factors can only be calculated 
numerically and the following discussion about the finite size behavior of SfP'(q,cu) has 
to be based on these numerical data. In figure [22] we show the sum of all two-particle 
form factors squared, denoted as l2{q), as a function of inverse length for different fixed 
momenta q and chain lengths up to N = 2600. On this scale it is hardly possible to 
detect any finite size effects at all. If we plot each curve for a fixed q- value separately on a 
finer scale, however, we see that the two-particle contribution is decreasing slightly with 
increasing length (see figure [23!) . By using extended precision arithmetics we checked 
that this slow decrease of the two-particle contribution with increasing length is not 
a numerical artefact. Based on the fits shown in figure [23] it is impossible to decide 
whether the two-particle contribution decreases (possibly logarithmically) to a finite 
value or vanishes with a power law. If it does vanish with a power law, the exponent is 
extremely small and apparently also q dependent. 

Although we do not have a detailed explanation for the observed size dependence, 
the following scenario seems plausible: Throughout the gapless regime, two-particle 
(defined in the Bethe Ansatz sense of shifted quantum numbers; these do not correspond 
to two-particle states in the field theory) contributions vanish in the thermodynamic 
limit at finite field, but are finite at zero magnetic field. At finite field (in contrast to zero 
field), there is room in the choice of quantum number configurations to accommodate 
multiparticle excitations yielding solutions to Bethe equations in terms of real rapidities. 
Within the gapless regime, backflows produced by finite numbers of excitations at the 
Fermi boundary vanish. With the possibility of adding particles without modifying 
the state substantially, it is probable that the contributions from 2, 4 and any finite 
number of excitations vanish in the thermodynamic limit. We can expect to get a good 
approximation for the exact thermodynamic result by summing over finite numbers of 

| Here, we call 'particle' an excitation obtained in the Bethe Ansatz by adding or removing a quantum 
number from the ground state. For the sake of simplicity, we treat bound states in a simplified manner 
and include their contribution in the appropriate term (i.e.: a two-string bound state counts as two 
excitations). 
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Figure 22. Sum of F 2 (q,u>) for all single particle-hole as a function of inverse length 
for different momenta q = 27r/20, • • • , 27r/50. Here A = 0.25 and a = —0.1. 
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Figure 23. The same as figure [22] The blue curves are fits ~ a + b/\n(l/L) 
with fit parameters a,b. The fit parameters (for decreasing q) are given by a = 
0.04321,0.034687,0.021739,0.0174 and b = -0.00245,-0.00111,-0.00023,-0.00011. 



The red curves are fits 



with fit parameters a = 0.0439,0.035,0.0218,0.0174 



and b = 0.00108, 0.00061, 0.00020, 0.00012. 
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excitations for a finite chain, but for very large values of N more and more families 
of states have to be included in the sum. In other words, although the correlation 
weight distribution among families of excited states (which is not an observable) might 
show some highly nontrivial size dependence, the full correlation function (the true 
observable here), which is obtained by summing over all states, shows a much weaker 
size dependence. This is supported by the analysis in (20] of the elementary blocks of the 
correlation functions for the finite and infinite chains. In the gapless regime all blocks 
(from which any correlation function can be obtained) differ from their thermodynamic 
limits by terms of order 1/N. 

This scenario is corroborated by the agreement found when comparing the 
numerical Bethe ansatz results to the field theory formulas for the high-frequency tail 
as well as by the evaluation of the sum rules in the previous subsection. The Bethe 
ansatz data in figure [19] correspond to four-particle contributions for a chain of length 
N = 600. 6, 8, ... and higher particle contributions have been neglected. So the four- 
particle contribution for a chain of length N = 600 is apparently a good approximation 
to the thermodynamic limit result. If each n-particle contribution does indeed vanish in 
the thermodynamic limit this suggests that it happens in such way that «Sjv(?,w) does 
not change significantly in the tail region. We call this the "compensation scenario" 
because when Sfj- starts to go to zero, S% comes in and "compensates" Sff. Similar 
compensations occur with S%, S$, etc. 

Concerning the peak region, we also believe that Sf? as shown in figure [T3l is a good 
approximation of the lineshape in the thermodynamic limit. It nearly saturates the sum 
rules indicating that even for N = 6000 the contribution of higher particle states are 
negligible. Furthermore, the width of the peak obtained from Sfj in section 15.21 does 
agree with the field theory result in section 15.11 If Sfj indeed vanishes for iV — > oo it 
again seems to get compensated by higher particle excitations in such a way that the 
lineshape does not change significantly. A notable exception might occur very close to 
where the methods of [28] predict a very singular form for S(q,u) caused by 
large numbers of excitations very close to the Fermi energy. It is difficult to determine 
whether or not S%(q,u) approaches this form at large N and it seems to be possible 
that SN(q,uj) could continue to change near ul,u(q), becoming more singular, out to 
very large N. 

9. Conclusion 

Based on a low-energy effective theory which includes the leading (band curvature 
type) irrelevant operators we have studied the longitudinal dynamical structure factor 
S zz (q,uo) for the XXZ spin- 1/2 chain in a magnetic field. By comparing with results 
for free fermions we have conjectured a method to sum up the entire perturbation series 
in one of these irrelevant operators allowing us to obtain an approximation for the 
shape of the peak of S zz (q,u) which is valid for small q and is non-perturbative in the 
interaction strength (anisotropy A). Besides the velocity v and Luttinger parameter K, 
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the important parameters to determine the lineshape are the coupling constants r/± of 
the leading (dimension-three) irrelevant operators, which we relate to derivatives of v 
and K with respect to the magnetic field. A summation of the entire series is necessary 
because perturbation theory in the band curvature terms is divergent on shell, uj ~ vq, 
although these operators are formally irrelevant. Our field theory approach is valid in 
the regime jq 2 <C rj_q <C v, where 7 is of order of the coupling constants of the next- 
leading (dimension-four) irrelevant operators which we neglected for the case of a finite 
magnetic field. The result is a box shaped peak with width 5u q = \rj-\q 2 and height 
K/\r)_\q, similar to the exact S zz (q,u) for the XX model (free fermion point). 

The field-theoretical results for the width of the peak are supported by Bethe Ansatz 
calculations. Since the XXZ model is integrable, we used the Bethe Ansatz equations 
in the thermodynamic limit to determine the parameters r)±(A,a) numerically, so that 
the low-energy effective theory and the obtained results for S zz (q,u) are parameter 
free. We have shown that the width of the peak obtained in field theory agrees with 
the analytically calculated width of the two-particle continuum in the Bethe Ansatz. 
Furthermore, we have demonstrated that the form factors obtained numerically by Bethe 
Ansatz approach the flat distribution predicted by field theory for q — > 0. Applying our 
results to the strongly interacting case (large A), we found that for A > cos(7r/5) ~ 0.81 
the parameter r]_ goes through zero for a finite value of the magnetic field. At the 
"inversion point" where r)-(A,ai nv ) = the q 2 scaling breaks down and the width 
of the two-particle continuum scales like 5uj q ~ q 3 . The q 2 scaling is recovered for 
< \a\ < o~i riv . As a result, the width Su q is a non-monotonic function of a, with a 
minimum at the inversion point. 

The power-law singularities found in [28] for A <C 1 near the lower and upper 
thresholds u>l{q), ^u{q) are 110 1 captured by our calculations. Within the effective low- 
energy theory these singularities seem to be related to higher dimension operators. 
Within the Bethe Ansatz, on the other hand, it is not clear if these singularities can 
be obtained by considering only the form factors for two-particle states. It seems 
possible that the finite size effects near these boundaries are complicated and form 
factors for multi-particle excitations at very large system sizes have to be studied. 
Nevertheless, for finite chains and small A the behavior of the dominant form factors 
for the two-particles states agrees qualitatively with the result of [28]. However, taking 
into account the energy dependence of the density of states leads to a maximum and 
a minimum of S zz (q,u) inside the two-particle continuum. In the strongly interacting 
regime A > cos(7r/5), we have found that for |er| < a inv the dynamical structure factor 
exhibits a rather distinct lineshape, reminiscent of a converging X-ray singularity at the 
lower threshold. 

We also showed that in the interacting case S zz (q, uj) has a high-frequency tail 
5S zz (q, uj). Within the effective theory this tail is related to an irrelevant operator 
(with coupling constant 77+) mixing excitations at the right and left Fermi points. 
Contrary to the calculation for the on-shell region, this term can be treated in finite- 
order perturbation theory for 5u q < w - < J and we find that the tail for finite 
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field decays as SS zz (q, u) ~ q A / {u 2 — v 2 q 2 ). This result is again supported by numerical 
calculations based on the Bethe ansatz. 

We have proposed that the integrability of the XXZ model is manifested in the low- 
energy effective Hamiltonian at the order of the dimension-four, band curvature type 
operators. The conservation of the energy current operator imposes that the interaction 
denoted as ( 3 is absent. This has consequences for the lineshape of S zz (q,u) at zero 
magnetic field, since a nonzero £3 would change the behavior of the tail near the upper 
threshold of the two-particle continuum. 

One promising test of our theory would be to measure, by means of inelastic neutron 
scattering experiments, the width of the peak as a function of q and h (equations (15.141) 
and (I3.69P ). In some spin-1/2 compounds it is experimentally possible to go all the way 
up to the saturation field. The main limitation is the low intensity of the signal for 
small-g scattering. However, one important point that may facilitate the experiment 
is that in the transverse channel the low energy spectral weight is shifted to a finite 
wave- vector, ±<r (the magnetization). The gap at q = is of order h for transverse 
excitations. So there is a "protected region" of zero transverse spectral weight at small 
q and u> inside of which the longitudinal structure function could perhaps be observed. 

Finally, we would like to emphasize that the formulas (I3.69P and (13.701) for the 
coupling constants of the irrelevant operators are also valid for nonintegrable models. 
This allows us to predict the width of the dynamical structure factor S zz (q,u) at small 
q once the field dependence of v and K is determined from thermodynamic quantities. 
For example, the value of (Ti nv , below which we expect to see nontrivial effects due to 
strong interactions, can be increased by adding a ferromagnetic next-nearest neighbour 
interaction. The isotropic J\ — J2 model also contains a marginally irrelevant operator, 
whose amplitude can be tuned by the Ji interaction. A ferromagnetic J2 (J2 < 0) 
would increase the constant do inside the logarithm in equation (I5.8ip . On more general 
grounds, the nonmonotonic behavior of the width 5u q and the inversion of the lineshape 
should occur whenever the derivatives of the velocity and the Luttinger parameter with 
respect to magnetic field/chemical potential have opposite signs and the latter one is 
singular. In principle, this could also be observed in the dynamical structure factor of 
quantum wires, since OK/ dn, where n is the electron density, changes sign and diverges 
in the low-density limit (Wigner crystal regime) |62j . This suggests that the evolution 
of the lineshape as a function of density could be richer than what was proposed in [28J . 

Open questions that we leave for future work include the behavior of S zz (q, u) near 
the upper and lower threshold, which will require a careful analysis of the interplay 
between more irrelevant operators in the effective model, a field theory calculation of 
the width for zero field, the thermodynamic limit of the form factors calculated in the 
Bethe Ansatz and further effects of integrability on the lineshape. 
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Appendix A. High-frequency tail for the zero field case 
In this appendix we derive the results ( 17.23D and ( 17.24p . 



Tail from £ + interaction 



Since the ( + vertex in (17.71) has two R and two L legs, the correction to x{Qi ^0 is 
separable into 8\ = &Xrr + $Xll, where 



K 



In L 



Il RLL (q,iu). 



(A.1) 



n_RLL is the bubble with one right- and two left-moving bosons (first diagram of figure 
IA1I) given by 



U RLL (q,iuj) 



dxe 



iqx 



dTe iu)T D ( °\x,T) 



D$\q,r) .(A.2) 



The expression for 8\ll is obtained from (lA.ip and (1A.2[) by exchanging R <-> L. After 
doing the Fourier transform and integrating over the internal frequencies, we find 

2ti 2 CI ^ k 1 k 2 {q + k l + k 2 ) 



U RLL (q,iuj) 



E 



L 2 ' iu — vq — 2v(k\ + k 2 ) 
fei,fe 2 >o y v 1 



(A.3) 



where k\ t2 — 2irni i2 /L, with n\ j2 integers. Taking the imaginary part of the retarded 
self-energy, we have 



Im U& L (q,u>) 



2nXl f2n 
vL 2 



x nin 2 {n + n\ + n 2 )5(£ — n — 2ni — 2ra 2 )(A-4) 

ni,ri2>0 

where we have used q = q n = 2irn/L and uj = uj£ = 2-kvE/L. Notice that this implies 
that the energy levels in the tail are discrete and separated by Anv/L. We evaluate the 
sum on the righthand side of (1A.4|) as follows 

nin 2 {n + n\ + n 2 ) 5(£ — n — 2n± — 2n 2 ) 



ni,n2>0 
oo 



y^(n + m) 5(£ — n — 2m) ni(m — n\) 



m=l 

oo 



ni=0 



Em 3 (n + m) / 1 \ 
6 1 ~ m? J 



— n — 2m) 
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Figure Al. Diagrams at 0(C+) for the calculation of the tail. 



6 4^ 



n\ ( £ + n 



£-n 



2irv 



5(UJ - LUi), 



(A.5) 



with Ui = 2ttv£/L, £ = n + 2,n + A, ... . Substituting (IA.5D in (IA.4[) and using flA.ip . 
we find 



-2Imx^(g,^ 



192w 2 V £ 

^n 2 (£ 2 -n 2 ) 



x 



6{u-u t ). (A.6) 



Likewise, we have 



2ImxS(<z,a, 



192v 2 V £ 



X 



1 - 



6(u-u t ). (A.7) 



Finally, the contribution of the ( + interaction to the high-frequency tail is 
6Sg(q,u) = ^^F 2 + (g n ,^)<5(^-^), 



(A.8) 



with 



F 2 + (q n ,u e ) 



K(l /27T 



n 2 (£ 2 -n 2 ) 



96w 2 V ^ 

In the thermodynamic limit L — > oo (and £ ± n ^> 1), we obtain 

^(C + H 2 2 (tf - vV 



192v 



(£-n) 2 (£ + n) 2 _ 
)btaii 

6{uj — vq). 



. (A.9) 



(A.10) 



Tai/ from Ai interaction 

The perturbation theory in the Umklapp interaction for a finite system requires that we 
treat the zero mode operators. So consider the expressions for spin operators 

It { , \ 

(A.11) 
(A.12) 



Sz 
j ~ \/ — "X 



+ (— l) j const x cos ( V AnK(j) 
S7 ~ const x e -*V*7* fl 



-1) J + cos (v47rK(J) 



Periodic boundary conditions for the spin operators imply that we can regard and 9 
as compactified fields with radius R = (1/AttK) 1 ^ 2 and R = (K/tt) 1 ^ 2 , respectively. In 
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(A.13) 
(A.14) 



where S z and m are integers. In a finite system with periodic boundary conditions we 
use the mode expansion for the bosonic fields 



vt x 

4> (x, t) = 4> + n — + Qo- 



Rg—iqnivt—x) 



+ a^e- iq " {vt+x) + h.c] , (A. 15) 



n>0 



where q n = 2irn/L. The operators <fro and Ho are associated with the zero mode and 
satisfy [4>q, ILq] = i. The compactification of <p quantizes the eigenvalues of Qq to be 
S z y/n/K. It follows from Eqs. flA.111) and (1A.15I) that S z corresponds to the total spin 
in the chain. We shall be restricted to the subspace S z = 0, to which the ground state 
for even L belongs. From dt<fi = vd x 9, we get 



Rg—iqnfyt—x) 



+ a^e- iqn{vt+x) + h.c] , (A. 16) 



n>0 



with [#o>Qo] = i- The eigenvalues of Ho are then my/kirK, m integer. Therefore, for 
the Hamiltonian 



H = l -Jdx [(djf + (d x <pf] 



2L 

we obtain the spectrum (S z = 0) 
2ttv 



n>0 



E 



L 



? K + J2n( m K + m L n ) 

n>0 

where m^ ,L = 0, 1, 2, . . .. The corresponding wave function is 
|*) = exp UmV^K^o} ]J (a?)"* (a?)"* |0> . 



(A.17) 
(A.18) 

(A.19) 
(A.20) 



n>0 



Since translation by one site takes <p — > (p+irR and |\&) — > (— l) m \^f) [63J, this symmetry 
implies that only intermediate states with even m couple to the ground state via S z . 
For the Umklapp interaction defined in f)7.7j) . the O(Xf) correction to x(g, zu;) is 



5y(q, iu) = ( — 

A w ' ; 87T V27T 



dxe~ 



iqx 



dr e 1 ' 



d 2 X\ / d 2 X2 



x ^(xy^^^e-'^^d^)) + (1 <-> 2). (A.21) 
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(x) (0) e *V^(i) e -iV^(2)\ 



con 



= 16nK [(0(x)0(l)> (0(0)0(1)) - (0(x)0(l)> (0(0)0(2)) 

x < / e i4v / ^^(l) e -i4v / ^</.(2)\ 

As a result, 5\ can be cast in the form 



2)] 



/ 





2 


~D(°\q,iujy 


V 27T , 


) 





[n( g ,^)-n(o,o)] 



(A.22) 



(A.23) 



where 



U{q,iu) = ~J dxe~ iqx J dre™ (e iiV ^^ x ^ e -«V^(o,o)^ . ( A .24) 

The correlation function n(x,r) = ^ e * 4v/ ^^ x ' T )e _44v/ ^'^ ' " ) ^ for a finite system 
has to be calculated using the mode expansion (lA.15j) including the zero mode. Note 
that the operators in ( 1A.21I) couple the ground state to states with m = ±2, since we 
are calculating matrix elements of the form 

J2V4nK4> - 



d x <pe 



. 



(A.25) 



U(z,z) 



4K /-,\ < 



Since e l4 v / ^<^(z,2) j s a primary field of holomorphic weight (2K,2K), the correlation 
function in the infinite complex plane is given by 

T 

We use the "CFT normalization condition" of [H]. The correlation function for a finite 
system is obtained using the conformal mapping z = e 27T ^^ L , z = e 2n ^ L , where £ = vr+ix 
and £ = vt — ix, with < x < L. The result is 



(A.26) 



n(£,0 



n/L 



AK 



n/L 



sin7r(a; + ivr)/L 



AK 



. (A.27) 



sin7r(x — ivr)/L 

In order to calculate Im H ret at zero temperature, we switch back to real time with the 
prescription ivr —> vt — ia, a — > + . We then calculate 

f L r + °° T 7r / T 

U(q,u) = - i 1 J ~--- i<lx 1 





r+oo 


1 dxe~ iqx j 


1 dte iuJt 




— oo 



sin n(x — vt + ia) / L 



n/L 



-i AK 



X 7 " —TTT > ( A - 28 ) 

sin tt{x + vt — ia) I L 

which has the property IT(g, u) = 2TmIT ret (g, uj). We also use the fact that for a periodic 
function with discrete modes q n = 2nn/L 



(A.29) 
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5> 5 l27~ 



2n 



n . (A.30) 



So we will consider the integral on the entire plane and eventually cancel a sum over 
delta functions for the discrete momenta. Performing a change of variables 

'qL 



2vr 



n 



l /7T\ 



sk r+°° 



2v \LJ 

r+oo 

x / dx_e- i{vq+uj)x - /2v 



7i (u — ia) 

sin ■ 



dx,e i{uJ - vq)x+/2v 



TT I X-L- — ia] 



-I -AK 



sin 



7r (x_ + ia) 
sin - 



L 



L 

-AK 



(A.31) 



where x± = x ± vt, we are left with integrals of the form 

hoo r / M-4if 

h = I due iru 



+oo 



0{r) £ 



L 

dz e irz Disc 



BC n 



7T (z — ia) 



sin 



-AK 



(A.32) 



n=— oo 

where SC n is the branch cut z = nL + ia + iy, < y < oo and Disc f(z) = 
/(z — 0~) — /(z + + ) is the discontinuity of the function across the branch cut. By 
shifting to z' = z — nL, we get 

-AK 



= -e (r) eirnLe ~ 



dz e irz Disc 

BCo 



TT (Z — ia) 



sin 



We have 



Then 



Disc 



sinh 



Ti (z — ia) 



-AK 


sinh — 




Lj 



-AK 



2i sin 4ttK. 



. (A.33) 



(A.34) 



9 (r) 2i sminKe 



iAirK 



in(rL— AttK) 



x / dye~ rv 



smh — 



-AK 



We use 



t in{rL— AttK) 



£* 



rL - AttK 
2~k 



m 



and (from [26] 



ds [sinh (vrTs)] e 



nAK-l 
AK isz L 



B 2K - i , 1 - AK 



(A.35) 



(A.36) 



(A.37) 
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where B(x,y) = T(x)T(y)/T(x + y) is the Euler Beta function, and finally get 

2 4K L ( rL 

h = sin(4vrAn 9 (r) B \2K + — , 1 - 4 

7r \ 2n 



x 



E a 



rL - Att K 
2^ 



m 



Likewise, the integral 



du e 



sin 



7r (u + ia) 



-4K 



(A.38) 



(A.39) 



is given by 



2 4K L ( rL 

J sin(4vrA) 9 (f) B [2K + — , 1 - AK 
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2vr 



x 



rL — AttK 
2t\ 



m 



(A.40) 



In our case, r = (u — vq) /2v and r = (uj + vq) /2v. The two delta functions can be 
recombined to replace the second condition by q = q n = 2imj L. 



in 

= E* 



rL — 47rA 



— m 



2tt 



E* 



rL - AttK 



2ix 



Attv 



— m — 



2 *)EMg 



qL 



m 



(A.41) 



We can then cancel the second delta function and write 

^ 8K-2 

U(q,uj) = - Ai 



2vr 

T 



x B 2K 



sm 2 (AnK)B [ 2K 
(ui + vq) L 



(uj — vq) L 
Attv 



l-AK 



Attv 



/ 



where ug = 2nv(£ + AK) / L, £ — n, n+2, . . .. Finally, the contribution from the Umklapp 
operator to the high-frequency tail is 



2vr 



8S£(q,u) = —^2Ff 1 (q n ,wt)5(u-w t ), 



(A.43) 



with 



F 2 1 (q n ,u e ) = 2 



2\ l K\ 2 (2^ m ~ A 



7TV 



B AK 



L 

i-n 



— sin^(47rA) 



n 



1-AK)B[AK + 



(I 2 - n 2 Y 
+ n 



l-AK) . (AAA) 



In the limit L — > oo, we can use B(x, y) ~ r(y)x y for x — > oo and we find 

2X 2 K 2 



T 2 (AK) 



(2vf~ m q 2 (uj 2 - v 2 q 2 ) iK - 3 9(u - vq). 



(A.45) 
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Figure A2. Diagrams for the high-frequency tail involving the £3 interaction. 
Infrared- divergent tail for £3 ^ 

The tail of order vanishes at uj — > vq because the to dependence of Hrll (equation 
IA. 5|) cancels the factor of (uj — vq)~ 2 from the external legs of the corresponding diagram 
in Sxrr- It is easy to see that if the internal bubble had a different combination of R 
and L bosons (e.g., if we replaced Hrll by Hrrl) this cancellation would not happen, 
leading to a divergence at uj — > vq. The dangerous combinations of external legs and 
three-boson bubbles are excluded for the XXZ model, but are allowed for non-integrable 
models with £ 3 7^ 0. The extra diagrams that contribute to the tail (8uj q <C uj — vq <C J) 
to second order in the coupling constants are illustrated in figure IA21 The calculation 
is similar to the 0(£+) diagram. The result in the thermodynamic limit is 

" CI uj a + 6v 2 q 2 uj 2 + v i q 4 2(z( + u 2 + v 2 q 2 ~ 



v 2 v 2 (uj 2 — v 2 q 2 ) v 2 v 2 



.(A.46) 



Note the divergence as uj — > vq in the 0(Q) term. As in the finite field case, this 

divergence stems from the breakdown of perturbation theory in the band curvature 

terms at uj ~ vq. We expect that this £3 contribution, which again is only finite for 

nonintegrable models, smoothes out the behavior of S zz (q,uj) near the upper threshold 

Uu(q), where the high-frequency tail joins the on-shell peak. 
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